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Abstract. We discuss the algebro-geometric initial value problem for the 
Ablowitz-Ladik hierarchy with complex-valued initial data and prove unique 
solvability globally in time for a set of initial (Dirichlet divisor) data of full 
measure. To this effect we develop a new algorithm for constructing stationary 
complex- valued algebro-geometric solutions of the Ablowitz— Ladik hierarchy, 
which is of independent interest as it solves the inverse algebro-geometric spec- 
tral problem for general (non-unitary) Ablowitz-Ladik Lax operators, start- 
ing from a suitably chosen set of initial divisors of full measure. Combined 
with an appropriate first-order system of differential equations with respect 
to time (a substitute for the well-known Dubrovin-type equations), this yields 
the construction of global algebro-geometric solutions of the time-dependent 
Ablowitz-Ladik hierarchy. 

The treatment of general (non-unitary) Lax operators associated with gen- 
eral coefficients for the Ablowitz-Ladik hierarchy poses a variety of difficulties 
that, to the best of our knowledge, are successfully overcome here for the first 
time. Our approach is not confined to the Ablowitz-Ladik hierarchy but ap- 
plies generally to (1 + l)-dimensional completely integrable soliton equations 
of differential-difference type. 



1. Introduction 

The principal aim of this paper is an exphcit construction of unique global so- 
lutions of the algebro-geometric initial value problem for the Ablowitz-Ladik hi- 
erarchy for a general class of initial data. However, to put this circle of ideas 
into a proper perspective, we first very briefly recall the origins of this subject: 
In the mid-seventies, Ablowitz and Ladik, in a series of papers [3]-[6] (see also 
[1], [21 Sect. 3.2.2], [Zl Ch. 3]), used inverse scattering methods to analyze certain 
integrable differential-difference systems. One of their integrable variants of such 
systems included a discretization of the celebrated AKNS-ZS system, the pair of 
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coupled nonlinear differential-difference equations, 

-iat - (f - af3){a^ + a+) + 2a = 0, 

(f.f) 

-z/3t + (f - a/3)(/3- + /3+) - 2/3 = 

with a = a{n,t), [3 = P{n,t), {n,t) e Z x M. Here we used the notation f^{n) = 
f(n ± f ), n e Z, for complex-valued sequences / = {/(n)}„gz- In particular, 
Ablowitz and Ladik [4] (see also [71 Ch. 3]) showed that in the focusing case, where 
/9 = —a, and in the defocusing case, where fi = a, (jl.ip yields the discrete analog 
of the nonlinear Schrodinger equation 

-iat-(l±|ap)(a"+a+) + 2a = 0. (1.2) 

Since then there has been an enormous activity in this area and we refer, for in- 
stance, to [3 Ch. 3], [57], [S], [31], [33], [SSI, [S], [ID], 03], [13, [IS], [SD] and the 

extensive literature cited therein, for developments leading up to current research 
in this particular area of completely integrable differential-difference systems. Par- 
ticularly relevant to this paper are algebro-geometric (and periodic) solutions of 
the AL system (jl.ip and its associated hierarchy of integrable equations. The first 
systematic and detailed treatment of algebro-geometric solutions of the AL system 
(jl.ip was performed by Miller, Ercolani, Krichever, and Levermore [40] (see also 
[9], [To], [H], [IS], [39], [52]). Algebro-geometric solutions of the AL hierarchy were 
discussed in great detail in [32] (see also [25], [26], [53]). The initial value problem 
for the half-infinite discrete linear Schrodinger equation and the Schur flow were 
discussed by Common [T7] (see also [TS]) using a continued fraction approach. The 
corresponding nonabelian cases on a finite interval were studied by Gekhtman ^24] . 
In addition to these developments within integrable systems and their applications 
to fields such as nonlinear optics, the study of AL systems recently gained consider- 
able momentum due to its connections with the theory of orthogonal polynomials. 
Especially, the particular defocusing case /3 = a and the associated CMV matrices 
and orthogonal polynomials on the unit circle attracted great interest. In this con- 
text we refer the interested reader to the two-volume treatise by Simon |47| (see 
also [IH]) and the survey by Deift [TO] and the detailed reference lists therein. 

Returning to the principal scope of this paper, we intend to describe a solution 
of the following problem: Given p = (p_,p+) S Nq \ {(0,0)}, r € Ng, assume 
Q,(o) ^ ^(0) i^g solutions of the pth stationary Ablowitz-Ladik system s-ALp(a, /?) — 
associated with a prescribed hyperelliptic curve K-p of genus p = p- + p+ — 1 
(with nonsingular affine part). We want to construct a unique global solution a = 
a{tr), P = PiU) of the rth AL flow AL^(a, /?) = with a(io,r) = a^°^ , P{tQ,r) = (3^"^ 
for some to,r G Thus, we seek the unique global solution of the initial value 
problem 

AL,(a, P) ^ 0, 

s-ALp(aW,/?(°)) =0, 

where a = a{n, tr), (3 = (3{n, tr) satisfy the conditions in (|2.2p . 

Given the particularly familiar case of real-valued algebro-geometric solutions of 
the Toda hierarchy (see, e.g., [16], [29l Sect. 1.3], [5T1 Sect. 8.3] and the extensive 
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literature cited therein), the actual solution of this algebro-geometric initial value 
problem, naively, might consist of the following two-step procedureQ 

(z) An algorithm that constructs admissible (cf. Section 2]) nonspecial divisors 
'^fiin) G Sym^ ICp for all n G Z, starting from a nonspecial initial Dirichlet divisor 
T^P-ljio) £ Sym^/Cp. "Trace formulas" of the type (|3.30p and (|3.3ip (the latter 
requires prior construction of the Neumann divisor Dp from the Dirichlet divisor 
though, cf. (|4.3p ) should then construct the stationary solutions a'^'^\ /3'^*'-' of 
s-ALp(a,/3) = 0. 

(m) The first-order Dubrovin-type system of differential equations (|5.43p . aug- 
mented by the initial divisor 'Djx(na,ta ,.) — T^fi{no) (cf- step (i)) together with the 
analogous "trace formulas" p.30p , (|3.3ip in the time-dependent context should then 
yield unique global solutions a = a{tr), f3 — (3{tr) of the rth AL flow ALr{a, (3) — 
satisfying a(to,r) = a^°\l3{to,r) = ■ 

However, this approach can be expected to work only if the Dirichlet divisors 
T^ii{n,tr) £ Sym^/Cp, fl{n,tr) — {fli{n,tr), . . . , flp{n,tr)), yield pairwise distinct 
Dirichlet eigenvalues iJj{n,tr), j = for fixed {n,tr) G Z x R, such that 

formula (|5.43p for the time-derivatives fJ.j,t^, j — 1, . . . ,p, is well-defined. Analo- 
gous considerations apply to the Neumann divisors 2?^ G Sym^/Cp. 

Unfortunately, this scenario of pairwise distinct Dirichlet eigenvalues is not real- 
istic and "collisions" between them can occur at certain values of (n, tr) G Z x R. 
Thus, the stationary algorithm in step (i) as well as the Dubrovin-type first-order 
system of differential equations (|5.43p in step (ii) above breaks down at such values 
of (n,tr). A priori, one has no control over such collisions, especially, it is not 
possible to identify initial conditions 2?^(no,tor) some (?T.o,to,r) G Z x M, which 
avoid collisions for all (n, tr) € x R. We solve this problem head on by explicitly 
permitting collisions in the stationary as well as the time-dependent context from 
the outset. In the stationary context, we introduce an appropriate algorithm al- 
luded to in step (i) by employing a general interpolation formalism (cf. Appendix 
IB|) for polynomials, going beyond the usual Lagrange interpolation formulas. In the 
time-dependent context we replace the first-order system of Dubrovin-type equa- 
tions (|5.43p . augmented with the initial divisor I?/i(„o,f(, by a different first-order 
system of differential equations (|6.15p . (|6.23p . and (|6.24p with initial conditions 
(|6.25p which focuses on symmetric functions of /ii(n, tr), ■ ■ ■ , fJ^pin, tr) rather than 
individual Dirichlet eigenvalues fij{n,tr), j — 1, . . . ,p. In this manner it will be 
shown that collisions of Dirichlet eigenvalues no longer pose a problem. 

In addition, there is an additional complication: In general, it cannot be guar- 
anteed that fij{n,tr) and h'j{n,tr), j = l,...,p, stay finite and nonzero for all 
(n, tr) G Z X R. We solve this particular problem in the stationary as well as 
the time-dependent case by properly restricting the initial Dirichlet and Neumann 
divisors I'/i(„o,to r)^'^o{no,to r) ^ Sym^ ICp to a dense set of full measure. 

Summing up, we offer a new algorithm to solve the inverse algebro-geometric 
spectral problem for general Ablowitz-Ladik Lax operators, starting from a prop- 
erly chosen dense set of initial divisors of full measure. Combined with an appro- 
priate first-order system of differential equations with respect to time (a substitute 



^We freely use the notation of divisors of degree p as introduced in Appendix |A] (see also the 
beginning of Section [S]!. 
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for the Dubrovin-type equations), this yields the construction of global algebro- 
geometric solutions of the time-dependent Ablowitz-Ladik hierarchy. 

We emphasize that the approach described in this paper is not limited to the 
Ablowitz-Ladik hierarchy but applies universally to constructing algebro-geometric 
solutions of (l + l)-dimensional integrable soliton equations. In particular, it applies 
to the Toda lattice hierarchy as discussed in [30] • Moreover, the principal idea of 
replacing Dubrovin-type equations by a first-order system of the type (|6.15p . (|6.23|1 . 
and (|6.24p is also relevant in the context of general (non-self-adjoint) Lax operators 
for the continuous models in (1 -I- l)-dimensions. In particular, the models studied 
in detail in [28 can be revisited from this point of view. (However, the fact that the 
set in ()6.67p is of measure zero relies on the fact that n varies in the countable set 
Z and hence is not applicable to continuous models in 1 -t- 1-dimensions.) We also 
note that while the periodic case with complex-valued a, (3 is of course included in 
our analysis, we throughout consider the more general algebro-geometric case (in 
which a, (3 need not be quasi-periodic). 

Finally we briefly describe the content of each section. Section [2] presents a quick 
summary of the basics of the Ablowitz-Ladik hierarchy, its recursive construc- 
tion. Lax pairs, and zero-curvature equations. The stationary algebro-geometric 
Ablowitz-Ladik hierarchy solutions, the underlying hyperelliptic curve, trace for- 
mulas, etc., are the subject of Section [3] A new algorithm solving the algebro- 
geometric inverse spectral problem for general Ablowitz-Ladik Lax operators is 
presented in Section (4] In Section [5] we briefly summarize the properties of algebro- 
geometric time-dependent solutions of the Ablowitz-Ladik hierarchy and formulate 
the algebro-geometric initial value problem. Uniqueness and existence of global 
solutions of the algebro-geometric initial value problem as well as their explicit 
construction are then presented in our final and principal Section [6l Appendix |^ 
reviews the basics of hyperelliptic Riemann surfaces of the Ablowitz-Ladik-type 
and sets the stage of much of the notation used in this paper. Finally, various in- 
terpolation formulas of fundamental importance to our stationary inverse spectral 
algorithm developed in Section|3]are summarized in Appendix [B] These appendices 
support our intention to make this paper reasonably self-contained. 

2. The Ablowitz-Ladik Hierarchy in a Nutshell 

We briefly review the recursive construction of the Ablowitz-Ladik hierarchy and 
zero-curvature equations following [21] and [35] . 

Throughout this section we suppose the following hypothesis. 

Hypothesis 2.1. In the stationary case we assume that a, (3 satisfy 

a,f3e C^, a{n)P{n) ^ {0, 1}, n e Z. (2.1) 

In the time- dependent case we assume that a, (3 satisfy 

a(-,i),/5(-,0 eC^ ieM, a(n, •),/5(n, •) eCi(R), neZ, 

a{n, t)l3{n, t) ^ {0, 1}, (n, t) G Z x R. 

Here denotes the set of complex-valued sequences indexed by Z. For a dis- 
cussion of assumptions (|2.ip and (|2.2p we refer to Remark [331 

We denote by the shift operators acting on complex-valued sequences / = 
{/(n)}„ez G according to 

(5±/)(n)=/(n±l), nGZ. (2.3) 
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Moreover, we will frequently use the notation 

/± = S^f, f e C^. (2.4) 
Next, we define sequences {fi,±}eeNo, {.9«.±}^eNo, and {/i^±}feNo recursively by 

30,+ = 5Co^+, /o,+ = -co,+q;+, ho,+ = co,+/3, (2.5) 

9^+1,+ - 97+1,+ = ^K,+ + Ph,+^ ^ £ f^o, (2.6) 

f7+i,+ = fe,+ - (^{91+1,+ + .97+1,+)' ^ ^ ^0, (2.7) 

/i£+i,+ = h^^ + I3{gi+i,+ + 9j+i,+), ^ e No, (2.8) 

and 

50,- = ^co /o,- = co,-a, /lo - = -Co -/?+, (2.9) 

5^+1,- - .97+1,- = "/i^,- + Z^/r-, ^ e No, (2.10) 

/£+i,- = + a(.g£+i,- + .97+1,-), ^ e No, (2.11) 

V+i,- = - /3(.9€+i,- + 57+1,-)' ^ e No. (2.12) 

Here co,± G C are given constants. For later use we also introduce 

/-i,± = /»-i,± = 0. (2.13) 

Remark 2.2. The sequences {/f^+j^gNo, {5f,+}feNo, and {/i^^+j^gNo can be com- 
puted recursively as follows: Assume that 9^,+^ and hi,+ are known. Equation 
(|2.6p is a first-order difference equation in 9e+i.+ that can be solved directly and 
yields a local lattice function that is determined up to a new constant denoted by 
ci+i^+ S C. Relations (|2.7p and (|2.8p then determine /£+i,+ and hi+i,+, etc. The 
sequences {/<;,- j^eNo, {5^,-}feNo, and {/i^,_}£gNo are determined similarly. 

Upon setting 

7 = l-a/3, (2.14) 

one explicitly obtains 
/o,+ - co,+(-a+), = co.+ (-7+a++ + (a+)2/3) +ci,+ (-a+), 

50, + = |co,+ , 5i,+ = Co, + (-«+/?) + |ci,+, 

/io,+ = co,+/3, = co,+ (7/3" - a+/3^) +ci,+/3, 

/ - 2a+N (2.15) 
Jo,- = Co, -a, /i,-=co,-(7a - a p^j + ci,_a, 

50,- = |co,-, 51,- = co,-(-a/3+) + |ci._, 

V- =co,-(-/3+), = Co,- (-7+/3+++a(/3+)') +ci,_ (-/?+), etc. 

Here {cf.ij^gN denote summation constants which naturally arise when solving the 
difference equations for g^,± in (|2.6p . (|2.10p . Subsequently, it will also be useful to 
work with the corresponding homogeneous coefficients fe,±, ge,±, and he,±, defined 
by the vanishing of all summation constants Ck,± for k = 1, . . . and choosing 

co,± = 1, 

/o,+ = -a+, /o,--a, /£,± = ^,±|co.±=i,c,.±=o,j=i,...,f, ieN, (2.16) 
50, ± = 5' 5<?,± = 5£,±|co,±=i,cj,±=o,i=i,...,^, ^ e N, (2.17) 
/io,+ =/3, ho,-^-(3+, /i£,± = /i£,±|co,±=i,c,.±=o,j=i,...,^, ^gN. (2.18) 
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By induction one infers that 

ill 
fe,± = '^ci^k,±fk,±i ge,± = ^Q-fe,±5fe,±, hi^± = '^c£^k^±hk^±. (2.19) 

fe=0 fc=0 fc=0 

Next we define the 2x2 zero-curvature matrices 

(2.20) 



zP 1 
and 

/G-{z) -F-{z)\ 

for appropriate Laurent polynomials Fp(z), Gp{z), Hp(z), and Kp{z) in the spectral 
parameter z e C\{0} to be determined shortly. By postulating the stationary zero- 
curvature relation, 

= UVp-Vp+U, 





(2.22) 


relations 




0, 


(2.23) 


0, 


(2.24) 


0, 


(2.25) 


0. 


(2.26) 



In order to make the connection between the zero-curvature formalism and the 
recursion relations (|2.5p - (l2.12|) . we now define Laurent polynomials Fp, Gp, Hp, 
and Kp bj0 

p- P+-1 
Fpi^) = E fp-i.-^~' + E fp^-i-i.+z', (2.27) 

1=1 1=0 

p- p+ 
Gp{z) = Y.9p~-i,-^~' + T.Sp+-^.+^'' (2.28) 

1=1 1=0 

P--1 p+ 



Hp{z)= hp_.i.e^^z-' + J2hp+-i,+z^, (2.29) 

1=0 1=1 

p- p+ 

^pi^) = ^9p--L-z~'' + ^gp+-i,+z'^ = Gp{z) + gp^- -9p+,+- (2.30) 

1=0 1=1 

The stationary zero-curvature relation (|2.22p . = UVp — VpU, is then equivalent 

to 

-a(5p+,+ + .9p_,-) + fp+-i..+ - fp--i- = 0' (2.31) 



^In this paper, a sum is interpreted as zero whenever the upper Hmit in the sum is strictly less 
than its lower limit. 
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Thus, varying p± e No, equations (|2.3ip and p.32p give rise to the stationary 
Ablowitz-Ladik (AL) hierarchy which we introduce as follows 

- ' V/^(5p+,+ +5p-.-) + V-1.+ - V-i.- / (2.33) 

Explicitly (recalling 7=1 — a(3 and taking p_ = for simplicity), 

C(o,o)a 

s-AL(i,i)(a,/3) = «- , . ■ «+A . = 0' 



s-AL(o,)(a,/3)=( HO, 



s-AL(2,2)(a,^) 



7(co-a + co,+a+) - C(i4)a 
7(co,+/3" + co-f3+) + C(i^i)/3 

/-7(co,+a++7+ + Co -a 7" - a(co,+a+/3" + co-a^f3+)\ 
' -/3(co,^(a-)2 + co,+ (a+)2)) 

7(co,-/3++7+ + co,+/3— 7" - /3(co.+a+/3- + co,-a-/3+) 
V -a(co,+ (r)'+co,-(/3+)2)) / 



/-7(ci,_a +ci,+a+) -C(2,2)a\ 

7(ci,+r + ci,-/3+) + C(2,2)/3 y ' ' 



(2.34) 



represent the first few equations of the stationary Ablowitz-Ladik hierarchy. Here 
we introduced 

Cp = (cp_,_ + Cp+,+)/2, P± e No. (2.35) 

By definition, the set of solutions of (|2.33p . with p ranging in Nq and ci^± G C, 
£ £ No, represents the class of algebro-geometric Ablowitz-Ladik solutions. 

In the following we will frequently assume that a, (3 satisfy the pth stationary 
Ablowitz-Ladik system supposing a particular choice of summation constants C£^± G 
C, ^ = 0, . . . ,p±, p± e No, has been made. 

In accordance with our notation introduced in (|2.16p - (|2.18l) . the corresponding 
homogeneous stationary Ablowitz-Ladik equations are defined by 

s-ALp(a,/3) =s-ALp(a,/3)|^^_^^^_^^^^^g ,^^^^^^^^^, £ = (p_,p+) e Ng. (2.36) 

One can show (cf. [311 Lemma 2.2]) that = 9p~,- up to a lattice constant 

which can be set equal to zero without loss of generality. Thus, we will henceforth 
assume that 

5p+,+ =5p-,-, (2-37) 

which in turn implies that 

Kp = Gp (2.38) 

and hence renders Vp in (|2.2ip traceless in the stationary context. (We note 
that equations (|2.37p and (|2.38p cease to be valid in the time-dependent context, 
though.) 

Next, still assuming (|2.ip and taking into account (|2.38p . one infers by taking 
determinants in the stationary zero-curvature equation (j2.22p that the quantity 

Rp^Gl~ FpHp (2.39) 

is a lattice constant, that is, Rp — R^ = 0. Hence, Rp{z) only depends on z, and 
assuming in addition to (|2.ip that 

co,± e C \ {0}, (2.40) 
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one may write Rp sJ^ 

2p+l 

Rp{^^) - {clJA)z''P- l[{z~ E^), {i?„}^;o' C C \ {0}, (2.41) 

m—O 

where 

p = p-+p+-l. (2.42) 
In addition, we note that the summation constants ci,±, . . . , Cpj_^± in (|2.33p can 
be expressed as symmetric functions in the zeros Eq,. . . ,i?2p+i of the associated 
Laurent polynomial Rp in (|2.4ip . In fact, one can prove (cf. [3T]) that 

Q,± = co,±q(^±1), ^ = 0,...,p±, (2.43) 

where 

c,(^±i) (2.44) 



k 



(2jo)! ■ ■ ■ (2j2p+i)! ±jo p±i2p+i 



^ ^ 22fe(jo!)2 . . . (j2p+i!)2(2jo - 1) • • • (2J2P+1 - 1^^° " ■ ■ ■^2^+1 



A: e N, 



jo,--- J2p+1=0 
joH |-j2p+l = fc 

are symmetric functions of E^^ — {E^^ , . . . , -E^^+i). 

Next we turn to the time-dependent Ablowitz-Ladik hierarchy. For that pur- 
pose the coefficients a and (3 are now considered as functions of both the lattice 
point and time. For each system in the hierarchy, that is, for each p G N§, we 
introduce a deformation (time) parameter tp e ]R in a,/3, replacing a{n), f3{n) by 
a{n,tp),l3{n,tp). Moreover, the definitions ([f!20| . and ([2:27| ~ ([230l) ofU,Vp, 

and Fp,Gp, Hp, Kp, respectively, still apply. Imposing the zero-curvature relation 

Ut^ + UVp-V+U = 0, peNl (2.45) 
then results in the equations 

at^ = i{zFp + a{Gp + Kp) - Fp), (2.46) 

Pt^ = -i{PiGp + Kp) - Hp + z-'Hj,), (2.47) 

= z{G- - Gp) + zpFp + aH-, (2.48) 

= + aHp + Kp- Kp. (2.49) 

Varying p e Ng , the collection of evolution equations 

- ' V-«/3tp+/3(5p+,+ +-9p-,-)- V-i,-+ (2.50) 

tpeR,p^ ip-,P+) e No, 

then defines the time-dependent Ablowitz-Ladik hierarchy. Explicitly, 



V~*ft(o.o) +C(0,0)P/ 



■^We use the convention that a product is to be interpreted equal to 1 whenever the upper Hmit 
of the product is strictly less than its lower limit. 
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V"*/^t(i,i) +7(co,+/3 +co -/?+) +C(i,i)/3y 
AL(2,2)(a,/3) (2.51) 

/-*"i(2 2) - 7(co,+a'*'^7"'" + Co. -a 7~ - a(co, + co.-a"/3+)\ 
-/3(co,-(a-)2 + co.+ (a+)2)) 

"^/3t(2 2) + 7(co -/3^^7^ + co.+/? 7" - /3(co,+a+/3" + cq _a"/3+) 
-a(co,+ (r)'+co,-(/3+)')) / 



-7(ci^_a + ci,+a+) - C(2,2)Q! 
7(ci,+/3-+ci,_/3+)+C(2,2)/3 



0, etc. 



represent the first few equations of the time-dependent Ablowitz-Ladik hierarchy. 
Here we recall the definition of Cp in (j2.35|l . 

The special case p = (1, 1), co,± = 1, and C(i^i) = —2, that is, 



-iati ~ li'^ + Q!^) + 2a 
-z/?t, +7(r +/3+)-2/3 



= 0, (2.52) 



represents the Ablowitz-Ladik system (|l.ip . 

By (|2.50p . (|2.6p . and (|2.10p . the time derivative of 7 = 1 — a/3 is given by 

7tp = n{igp+,+ ~9p+,+) - {gp-,- -3jL,-))- (2-53) 

The corresponding homogeneous equations are then defined by 

AL^(a,/3) = AMa,/3)|^^^^^^^^^^^g ,^^_^^^^^^ =0, p = (p-,p+) e N^. (2.54) 

Remark 2.3. From (|2.23p - (|2.26p and the explicit computations of the coefficients 
ft,±^ gi,±, and he^±, one concludes that the zero-curvature equation (|2.45p and 
hence the Ablowitz-Ladik hierarchy is invariant under the scaling transformation 

a^aa = {ca{n)}nez, P ^ = {P{n) / c}nez, c £ C \ {0}. (2.55) 

In particular, solutions a, /3 of the stationary and time-dependent AL equations are 
determined only up to a multiplicative constant. 

Remark 2.4. [i) The special choices (3 = ±a, co.± — 1 lead to the discrete 
nonlinear Schrodinger hierarchy. In particular, choosing C(i,i) = —2 yields the 
discrete nonlinear Schrodinger equation in its usual form (see, e.g., [7i Ch. 3] and 
the references cited therein), 

- iat - (It + "+) + 2a = 0, (2.56) 

as its first nonlinear element. The choice j3 — ais called the defocusing case, (3 = —a 
represents the focusing case of the discrete nonlinear Schrodinger hierarchy. 
(ii) The alternative choice /? = a, co,± = leads to the hierarchy of Schur flows. 
In particular, choosing C(i i) = yields 

at - (1 - |an(a+ - a") = (2.57) 

as the first nonlinear element of this hierarchy (cf . [11] , [22] , [23] , [34] , [41] , [49] ) . 

Finally, we briefly recall the Lax pair {L,Pp) for the Ablowitz-Ladik hierarchy 
and refer to [33] for detailed discussions of this topic. 
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In the standard basis {5m}mez in defined by 



1, m — n, 
0, m ^ 



(2.58) 



the underlying Lax difference expression L is given by 
/ 



L 







-q(0)p(-1) -/3(-l)a(0) -a(l)p(0) p(0)p(l) 

p(-l)p(0) /3(-l)p(0) -/3(0)a(l) /3(0)p(l) 

-a(2)p(l) -/3(l)a(2) -q(3)p(2) p(2)p(3) 

p(l)p(2) /3(l)p(2) -/3(2)a(3) /3(2)p(3) 



V 







(2.59) 



P P Seven S + iP P Seven - a^pSodd)S - /3a+ 

+ (/3p+ <5evc„ - a++p+ SeM)S+ + p+ p++ <5odd S++, 



(2.60) 



where 5cvcn and Sodd denote the characteristic functions of the even and odd inte- 
gers, 

Seven =X-2Z^ Sodd — ~ Seven — X-2Z+1- (2-61) 

In particular, terms of the form — /9(n)a(n+l) represent the diagonal (n, n)-entries, 
n S Z, in the infinite matrix (|2.59|) . In addition, we used the abbreviation 

p = 7i/2 = (1 _^^)i/2^ (2.62) 

Next, let T be a bounded operator in the Hilbert space ^^(Z) (with scalar product 
denoted by ( • , • )). Given the standard basis ()2.58p in ^^(Z), we represent T by 



r=(T(TO,n)) 2, T(m,n) = (J„,r5„), (m, n) G Z^ 



(2.63) 



Actually, for our purpose below, it is sufhcient that T is an A^-diagonal matrix for 
some TV € N. Moreover, we introduce the upper and lower triangular parts T± of 
Tby 



T± = {T±{m,n)) 



2, T±{m,n) 



T{in, n), ±(n — m) > 0, 
0, otherwise. 



(2.64) 



Then, the finite difference expression Pp is given by 



^~ 2 



_((i^)+-(Lr)-^E^p-^-((^ 



1=1 



- -^CpQd, p = ip-,p+) e No, 
with Qd denoting the doubly infinite diagonal matrix 

g. = ((-i)%,.),,,,. 



(2.65) 



(2.66) 



and Cp = (Cp-,- +Cp^_+)/2. The commutator relations [fp,i] = and Lt^ — 
[Pp,L] = are then equivalent to the stationary and time-dependent Ablowitz- 
Ladik equations (|2.33p and (|2.50p . respectively. 
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3. Properties of Stationary Algebro-Geometric Solutions 

OF THE AbLOWITZ-LADIK HIERARCHY 

In this section we present a quick review of properties of algebro-geometric so- 
lutions of the stationary Ablowitz-Ladik hierarchy. We refer to [31] and [32] for 
detailed presentations. 

We recall the hyperelliptic curve /Cp of genus p, where 

2p+l 

/Cp : Tp{z, y)=y^- ^c^^z'^p- Rp(z) = - l[{z ~ E„,) = 0, 

m=0 

2 2p+l 



(3.1) 



Throughout this section we make the assumption: 

Hypothesis 3.1. Suppose that 

a,P eC^ and a(n)/?(n) ^ {0, 1} for all n e Z. (3.2) 

In addition, assume that the ajfine part of the hyperelliptic curve /Cp in (j3.1[) is 
nonsingular, that is, suppose that 

Em Em' for m m' , m,m' = 0,1, ... ,2p + 1. (3-3) 

The curve JCp is compactified by joining two points Poo±, Poo+ ^ -Poo_, but 
for notational simplicity the compactification is also denoted by /Cp. Points P on 
l^p\{Poo+ , Poo^ } are represented as pairs P = (z, y), where y(-) is the meromorphic 
function on /Cp satisfying Tpiz, y) = 0. The complex structure on /Cp is then defined 
in the usual way, see Appendix [^ Hence, /Cp becomes a two-sheeted hyperelliptic 
Riemann surface of genus p in a standard manner. 

We also emphasize that by fixing the curve /Cp (i.e., by fixing Eq, . . . , i?2p+i), 
the summation constants ci^±, . . . , Cp_,_^± in fp±,±, gp±,±, and /ip±,± (and hence in 
the corresponding stationary s-ALp equations) are uniquely determined as is clear 
from (|2.44p which establishes the summation constants Ci^± as symmetric functions 
oi E^^,...,E2p+i. 

For notational simplicity we will usually tacitly assume that p G N. 

We denote by {Mi("')}j=i,...,p and {i^j («)}j"=i,...,p the zeros of ( • Fp{ ■ , n) and 
( • y-^^Hp{ • ,n), respectively. Thus, we may write 



Fp{z) = -co.+a+z-"- \[{z 11,), (3.4) 
p 

Hp{z) = c„^+Pz-P-+' Hiz - V,), (3.5) 



and we recall that (cf. p.39p ) 

Rp — Gp — —FpHp. (3-6) 

The next step is crucial; it permits us to "lift" the zeros fij and Vj from the complex 
plane C to the curve /Cp. From p.6p one infers that 

Rp{z) - Gp[zf- =0, z e [nj, yk}i,k=i,...,p- (3.7) 
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We now introduce {Aj}j=i,- -,p C KLp and {^'j}j=i,...,p C /Cp by 

jjj{n) ^ {fij(n),{2/cQ,+)fij{n)P-Gp{fij{n),n)), j = l,...,p, n e Z, (3.8) 

and 

i>j(n) = (i^j(n), -(2/co,+)^'j(n)^-Gp(^'j(n),n)), j = 1, . . . ,p, n £ Z. (3.9) 
We also introduce the points Po.± by 

^2 2p+l 

Po,± = (0, ±(co,-/co,+)) e /Cp, ^ = n (3.10) 

0,+ m=0 

We emphasize that Po_± and Poo± are not necessarily on the same sheet of ICp. 
Moreover, 

|±(co,_/co,+) + 0(C), P^Po.±, C^z. 

Next we introduce the fundamental meromorphic function on ICp by 

(co,+/2)z-P-y + Gp(z,n) 
= Fpi^) (^-^'^ 



-Hp{z,n) 



(3.13) 



(co,+/2)z-P-y-Gp(0,n)' 

-P = iz,y) e /Cp, n e Z, 

with divisor (</)( • , 71)) of </>(•, n) given by 

{<j}{-,n)) =2'po,-£(n) -^'p^.AW' (3-14) 
using (|3.4p and (|3.5p . Here we abbreviated 

A = {Ai, • ■ • {t'l, • ■ • ,!>p} e SymP(/Cp). (3.15) 

For brevity, and in close analogy to the Toda hierarchy, we will frequently refer to 
fi and as the Dirichlet and Neumann divisors, respectively. 

Given (/)(•, n), the meromorphic stationary Baker- Akhiezer vector ^'(•,n,no) 
on /Cp is then defined by 



«'(P,n,no) 



vV'2(P,n-,no)y 

[ n"'=„o+i + «("')r l^', n')) , n > no + 1, 
'(/'i(P, n,7io) = < 1, n^UQ, (3.16) 

in::L„+i (z+ «(n')r (i',n'))"\ « < no - 1, 

[ n:'=„„+i (^/3(n')r (^^,^^0-' + 1), n > rio + 1, 

'(/'2(P, "•,"0) = 0(-P,"o) S 1, n = nQ, 

[n:Ln+i {z(3{n')<t>-{P,n')-' + 1)"\ n < no - 1. 

(3.17) 

Basic properties of <j) and ^' are summarized in the following result. 
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Lemma 3.2 ([3T]). Suppose that a, (3 satisfy p.2|) and thepth stationary Ablowitz- 
Ladik system (j2.33p . Moreover, assume (|3.ip and (j3.3p anc? ^ef P — {z,y) e /Cj, \ 
{Poo_i_, Poo-, Pq,+ , Pq,-}- (n^no) G Z^. Then (f> satisfies the Riccati-type equation 



cj>{P)cj,- [P] - r (P) + z^{P) - zP, (3.18) 



Of' 

as well as 

HJz) 

m^{P*) = (3.19) 
^p(z) 

Gp{z) 

^{P)+cP{P*)^2-^, (3.20) 

4>{P) ^ cj,{P*) = c,.+z-P- ^ (3.21) 

The vector \l/ satisfies 

U{z)^-{P) = '^{P), (3.22) 

yp(z)*-(P) = ~{i/2)co.+z-P-y^-{P), (3.23) 

V'2(/',n,no) ='/'(^,^)^i(-P,«,no), (3.24) 

F (z n) 

MP,n,no)MP*,n,no) ^ z"'"" ' , r(n,no), (3.25) 

_ff (z n) 

V'2(P,n,no)^2(P*,n,no) = z"'"° 7 ' , r(n,no), (3.26) 

Fp{z, no) 

ipiiP, n, no)ip2{P*,n, no) + ipi{P*,n, no)'iJj2{P, n, no) (3.27) 

G„(z, n) 
-2z-"" y y {n,no), 
Fp[z,no) 

ipi{P,n,no)ip2{P*,n,no) - ipi{P* ,n,no)i>2iP,n,no) (3.28) 

= -co,+z"-"°-P- , r(n,no), 
Fp(z, no) 

where we used the abbreviation 





= "0 + 1 ' 


1 









n = rto, (3.29) 



Combining the Laurent polynomial recursion approach of Section [2] with p.4p 
and p.Sp readily yields trace formulas for fi_± and /i^ -t in terms of symmetric 
functions of the zeros and ly^ of {■)'''- Fp and {■)P-^^Hp, respectively. For 
simplicity we just record the simplest cases. 

Lemma 3.3 ([31j). Suppose that a, (3 satisfy p.2|) and thepth stationary Ablowitz- 
Ladik system (j2.33p . Then, 



-^ = (-l)^+^^n,,, (3.30) 
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=a+/3-7^-^. (3.33) 



^ = (-ir^^n-.' (3-31) 

f]^, 7+^-^, (3.32) 

p 

E 

i=i 

Remark 3.4. The trace formulas in Lemma 13.31 illustrate why we assumed the 
condition a{n)P{n) ^ for all n G N throughout this paper. Moreover, the fol- 
lowing section shows that this condition is intimately connected with admissible 
divisors Vp^jVp^ avoiding the exceptional points -Poo±,-Po,±- On the other hand, as 
is clear from the matrix representation p.59|) of the Lax difference expression L, 
if a(no)/3(?T-o) — 1 for some uq £ N, and hence p{no) = 0, the infinite matrix L 
splits into a direct sum of two half-line matrices L±{nQ) (in analogy to the familiar 
singular case of infinite Jacobi matrices aS^ + + & on Z with a(no) = 0). 

This explains why we assumed a{n)(3{n) ^ 1 for all n G N throughout this paper. 

Since nonspecial divisors and the linearization property of the Abel map when 
applied to and I?i> will play a fundamental role later on, we also recall the 
following facts. 

Lemma 3.5 ([31], [22]). Suppose that a,/3 satisfy (|3.2p and the pth stationary 
Ablowitz-Ladik system (|2.33p . Moreover, assume p.ip and p.3p and let n G Z. 
Let Vfi, ft — {/ii, . . . , fip}, and Vc^, y_ — {v_y, . . . be the pole and zero divisors 

of degree p, respectively, associated with a, P, and cj) defined according to p.Sp and 
(EH), that is, 

Aj(»^) = (A*i("), (2/co,+)Aij(n)P-Gp(Aij(n),n)), j = l,...,p, 

— [3 34} 

i>j(n) = (i^j(n),-(2/co.+)j/j(n)P-G'p(i^j(n),n)), j^l,...,p. 

Then and I?^(„) are nonspecial for all n G Z. Moreover, the Abel map 

linearizes the auxiliary divisors T>j^ and in the sense that 

^Qo(%«)) =aQo(%no)) + {n~nQ)Ap^_{Poo^), (3.35) 

^Qo(%«)) =aQo(%«o)) + ("-"o)i.Po,_(f'oo+), (3.36) 
where Qq G ICp is a given base point. In addition, 

^Qo (%(")) =^Qo(%«)) +4po,_(f'oo_). (3.37) 

For theta function representations of a and P we refer to p32] and the references 
cited therein. These theta function representations also show that 7(n) ^ {0,1} 
for all n G Z, that is, the second condition in p.2p is satisfied for the stationary 
algebro-geometric AL solutions discussed in this section provided the associated 
Dirichlet and Neumann divisors are admissible. 



4. The Stationary Algorithm 



The aim of this section is to derive an algorithm that enables one to construct 
algebro-geometric solutions for the stationary Ablowitz-Ladik hierarchy for general 
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initial data. Equivalently, we offer a solution of the inverse algebro-geometric spec- 
tral problem for general Lax operators L in (I2.60p . starting with initial divisors in 
general position. 

Up to the end of Section [3] the material was based on the assumption that 
a,/3 € satisfy the pth stationary AL system (|2.33p . Now we embark on the 
corresponding inverse problem consisting of constructing a solution of (|2.33p given 
certain initial data. More precisely, we seek to construct solutions a, /3 £ sat- 
isfying the pth stationary Ablowitz-Ladik system (|2.33p starting from a properly 
restricted set M.^ of admissible nonspecial Dirichlet divisor initial data 'Dfi(no) 
some fixed rig G Z, 

/t(no) = {/ii(no), • ■ • , Ap("o)} e A^o, Mo d Syu\^ [ICp) , 

~ (4.1) 

= (Atj(no), (2/co,+)Aij(rio)^"Gp(^j(no),no)), j = 

For convenience we will frequently use the phrase that a, (3 blow up in this man- 
uscript whenever one of the divisors P/i or Do hits one of the points Poo±) Po,±- 

Of course we would like to ensure that the sequences a, (3 obtained via our 
algorithm do not blow up. To investigate when this happens, we study the image 
of our divisors under the Abel map. A key ingredient in our analysis will be ()3.35p 
which yields a linear discrete dynamical system on the Jacobi variety J{ICp). In 
particular, we will be led to investigate solutions Vf^, Vo of the discrete initial value 
problem 

SOo(%«)) ==aQo(%no)) + ("-'^oUp^,._(-Poo+), ^^^^ 

A("o) = {Ai(?^o), • ■ . , Ap("o)} e Syu\^{Kp), 

respectively 

Hno) = {hino), . . . , t'p(no)} e SymP(/Cp), 

where Qo G JCp is a given base point. Eventually, we will be interested in solutions 
Vfi, of (|4.2p . (|4.3p with initial data I'/i(„o) satisfying (|4.ip and to be specified 
as in (the proof of) Lemma H?T1 

Before proceeding to develop the stationary Ablowitz-Ladik algorithm, we briefly 
analyze the dynamics of (|4.2p . 



Lemma 4.1. Let n G Z and suppose that 2?/i(n) is defined via ()4.2p for some divisor 
T^Kno) ^ Sym^'(/Cp). 

(i) If'Dfi(^n) is nonspecial and does not contain any of the points Pq.±, Poo±, md 
contains one of the points Po.±, Poc± , then contains Pq,- or Poo_ 

but not Poc+ or Po,+ - 

{ii) If T^p,(n) is nonspecial and 'Df^(^n+i) is special, then 2?^(n) contains at least one 
of the points Poo+; ^oo_ ^.n-d- one of the points Po.+; ^o.-- 

(iii) Item (i) holds if n + 1 is replaced by n — I, Poo+ by Poo- , and Po.+ by Po^_ . 

(iv) Items (i)-{iii) also hold for T>c,(^ny 

Proof, (i) Suppose one point in I?/i(„+i) equals Poo+ and denote the remaining ones 
by Then (g^l) implies that ag^ (X»^(„+i)) + Aqo{Poo+) = aQo(2'^(n)) + 

Apg _iPca+)- Since I^fi{n) is assumed to be nonspecial one concludes = 
^/i(n+i) + T-^Po ~^ contradicting our assumption on 'Df^(^,iy The statement for Po,+ 
follows similarly; here we choose Qo to be a branch point of ICp such that Aq^ (P* ) — 
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(ii) Next, we choose Qq to be a branch point of /Cp. If P^(„_(.i) is special, then it 
contains a pair of points {Q, Q*) whose contribution will cancel under the Abel map, 
that is, aQ^(P^(„+i)) = aQ^(I?^(„+i)) for some 2?^(„+i) e SymP^^ {ICp) . Invoking 
^ then shows that aQ„(P^(„)) = ag^ (%n+i)) + ) + ^Qo (^o,-)- As 

was assumed to be nonspecial, this shows that 2^/i(n) — 2'ji(n+i) + T^Pao + 
2?P(, _ , as claimed. 
(Hi) This is proved as in item (i). 

(if) Since X'i>(„) satisfies the same equation as 2?/i(„) in (14. 2|) (cf. (I3.36|) ). items 
also hold for □ 

We also note the following result: 

Lemma 4.2. Let n G Z anrf assume that a»^d 2'i>(n) are nonspecial. Then 

2^/i(n) contains Po,- if and only ifD^^n) contains Poo_ - Moreover, T^fi(n) contains 
Poo+ if and only ifDpi^^^ contains Po.+ . 

Proof. Suppose a point in I'/i(„) equals -Po,- a-nd denote the remaining ones by 

%„)• By (ESZD, 

fiQo = fiQo (%«)) + (^0,-) + Ap^ _ (Poo J = «Q„ (%„)) + Aq„ (Poo_ ). 

(4.4) 

Since ^^{n) is nonspecial, 'Dc.^n) contains Poo_ , and vice versa. The second state- 
ment follows similarly. □ 

Let us call the points Poo^, Pcxd_, Po,+ , Sind Po,- exceptional points. Then 
Lemma 14.11 yields the following behavior of 'Djx(n) assuming one starts with some 
nonspecial initial divisor I'^(„o) without exceptional points: As n increases, 
stays nonspecial as long as it does not include exceptional points. If an exceptional 
point appears, T^fi(n) is still nonspecial and contains Po,- or Poo_ at least once (but 
not Po,_|- and Poo+)- Further increasing n, all instances of Po,_ and Poo_ will be 
rendered into Po,+ and Poo+ , until we have again a nonspecial divisor that has the 
same number of Po,-|- and Poo+ as the first one had of Po,_ and Poo_ . Generically, 
one expects the subsequent divisor to be nonspecial without exceptional points 
again. 

Next we show that most initial divisors are well-behaved in the sense that their 
iterates stay away from Poo±, Po.±- Since we want to show that this set is of full 
measure, it will be convenient to identify Sym^(A^p) with the Jacobi variety J(/Cp) 
via the Abel map and take the Haar measure on J(/Cp). Of course, the Abel map 
is only injective when restricted to the set of nonspecial divisors, but these are the 
only ones we are interested in. 

Lemma 4.3. The set Aio C Sym^(/Cp) of initial divisors P^(„p) for which ^?/i(n) 
and 2?i>(„), defined via (j4.2p and (j4.3p . are admissible {i.e., do not contain the 
points Poo±, Po,±) o,nd hence are nonspecial for all rt S Z, forms a dense set of full 
measure in the set Sym^(/Cp) of positive divisors of degree p. 

Proof. Let A^oo.o be the set of divisors in Sym^(/Cp) for which (at least) one point 
is equal to Poo-t or Po,±. The image aQ^{M.aDfi) of A^oo,o is then contained in the 
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following set, 

aQ„(A^oo,o)C y UQ„(P)+«Q„(SymP-i(/C,))) C J(/Cp). (4.5) 

-Pe{Po,±,Poo±} 

Since the (complex) dimension of Sym^^^ (ICp) is p — 1, its image must be of measure 
zero by Sard's theorem (see, e.g., O Sect. 3.6]). Similarly, let Msp be the set of 
special divisors, then its image is given by 

ag„(A1sp) = aQ„ (SymP-2(/C,)), (4.6) 

assuming Qq to be a branch point. In particular, we conclude that aQ^^{A4sp) C 
Q:Q^{A4oc,o) and thus aQ^{M.sing) = aQo(-^oo.o) has measure zero, where 

A^sing =Xoo,0U TWep- (4.7) 

Hence, 

5^ - U {aQo (Msing) + nAp^ _ (Poo+ )) and 5. = + Ap^ _ {Poo_ ) (4.8) 

are of measure zero as well. But the set iS^ U S^, contains all initial divisors for 
which X'^(„) or X'i;(„) will hit Poo± or Pq.±, or become special at some n € Z. We 
denote by Aio the inverse image of the complement of the set 5^ U under the 
Abel map, 

Mo = aQj(Symf (/Cp) \ (5^ U 5,)). (4.9) 
Since Mq is of full measure, it is automatically dense in Sym''(A^p). □ 

Next, we describe the stationary Ablowitz-Ladik algorithm. Since this is a 
somewhat lengthy affair, we will break it up into several steps. 
The Stationary Ablowitz Ladik Algorithm: 

We prescribe the following data 
(i) The coefficient a(no) £ C \ {0} and the constant co,+ £ C \ {0}. 
(m) The set 

{Emf^+^ CC\{0}, Em ^ E^, for m ^ m\ m, m' = 0, . . . , 2p + 1, (4.10) 
for some fixed p G N. Given {E„i}m=0' introduce the function 

/ X 2 2p+l 

Rpi^)=(^) n(^-^'») (4.11) 

and the hyperelliptic curve ICp with nonsingular affine part as in p.ip . 
(in) The nonspecial divisor 

%„o) e SymP(/Cp), (4.12) 

where fl{nQ) is of the form 

Kno) = {Al("o),---,Al('^o),---,Ag(no)'---'A9(no)} (4-13) 

v ' " V ' 

Pi (no) times Pg(„g) (no) times 

with 

fj-kim) ^ {^ik{no),y{flkina))), ^ikim) ^ ^ik' {m) for k ^ k' , k, k' ^ 1, . . . , q{no), 

(4.14) 
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and 



9(no) 



Pfe(no) e N, k ^ 1,. . ., q{nn), ^ Pk{m) 



(4.15) 



fc=i 



With {Em}^t,ln 2?^(no)i ot{nQ)^ and co,+ prescribed, we next introduce the following 
quantities (for z G C \ {0}): 



, 2p+l \ 1/2 <j("n) 

a+(no) -a(no)( [] ^™ ) 11 A'fcK)-^''^""^ 

^ m=0 ^ fc=l 
2p+l 



W Era, 



m— 



Fp(z,7io) = -co^+a+(no)z-P- [] (z - A'fcK)^^"''^ 
1 / 1 



k=l 



Gp{z, no) 



2\a{no) a+{no)J - 

z 



Fp{z,no) 



(4.16) 

(4.17) 

(4.18) 
(4.19) 

P=(C,v)=iikino) 



(Pk(no)-i-l 
(l(^Pk(na)-l- 



k=l e=o 

lino) ^ 

fe'=l,fc'^fe / 



C=Mfc("o) 



Here the sign of the square root is chosen according to (I4.14p . 

Next we record a series of facts: 
(I) By construction (cf. Lemma FB.II . 



d'{Gp{z,noy 



dz^ 



d'^Rp(z) 



dzJ 



z=fj.k (no) 



z e C \ {0}, 1^0,... ,pk{no) - 1, fc = 1, . . . , q{no). 

(II) Since 'Djx{riQ) is nonspecial by hypothesis, one concludes that 

Pk{no) > 2 implies Rp{iik{no)) 7^ 0, fc = 1, . . . , q{no). 

(III) By (|419ll and one infers that Fp divides Gl - Rp. 

(IV) By (|i7TT|) and ((iTTg)l one verifies that 

Gp{z,nof - Rp{z) = 0(z'P+-i), 



Gp{z,nof - Rp{z) = 0(z-2''-+i 

— — 2^0 



)• 



By (III) and (IV) we may write 

Gp{z,nQf ~ Rp{z) = Fp[z,no)Hq^r{z,no), z e C \ {0}, 



(4.20) 

(4.21) 

(4.22) 
(4.23) 

(4.24) 



for some q G {0,...,p_ — 1}, r £ {0, where Hgr(z,no) is a Laurent 

polynomial of the form c^qZ^'^ + • • ■ + c^z''. If, in fact, i?o,o = 0, then Rp{z) = 
Gp{z, no)^ would yield double zeros of i?p, contradicting our basic hypothesis (|4.10p . 
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Thus we conclude that in the case r = q = 0, i?o.o cannot vanish identicaUy and 
hence we may break up (|4.24|) in the foUowing manner 

Gp{z,no) + ico.+/2)z-P-y ff„,.(z,no) 
(l>(P, no) = ^ - 



Fpiz, no) Gp{z, no) - {co,+/2)z-p- y ' (4.25) 

P = {z,y) e ICp. 

Next we decompose 

r+q 

Hg,riz,no)^Cz-'>l[{z-i^j{no)), ^eC\{0}, (4.26) 

where C G C \ {0} and {^'j(?^o)}jii C C (if r = g = we replace the product in 
(|4.26p by 1). By inspection of the local zeros and poles as well as the behavior near 
Po,±, Poo± of the function (/){ ■ , no) using (|3.1ip . its divisor, (0( • , no)) , is given by 

{(pi-,no)) ^'Dpo,^oino) -2^Poo_A(«o)' (4-27) 

where 

£l(no) = {Po Po-^hino), ■■■ ,Vr+q{no), Poo+, Poo+}- (4.28) 

^ . ' " . ' 

p_—l — g times p+—r times 

In the following we call a positive divisor of degree p admissible if it does not 
contain any of the points Poo± i ^o.±- 
Hence, 

2?i>(no) is admissible divisor if and only if r = p+ and q — — 1. (4.29) 
We note that 

«Qo(^^(no)) = «Qo(%«o)) +4po,J^ooJ, (4.30) 

in accordance with (|3.37p . 

(V) Assuming that (I4.22p . (|4.23p are precisely of order z^'-^p^^^\ that is, assuming 
r ^ pj^ and (7 = p_ — 1 in (|4.24p . we rewrite (|4.24p in the more appropriate manner 

Gp{z,no)^ -Rp{z) = Fp(z,no)Hp{z,no), ^eC\{0}. (4.31) 

(We will later discuss conditions which indeed guarantee that q — p^ — 1 and r = p+ , 
cf. (|4.29p and the discussion in step (X) below.) By construction, Hp{- ,no) is then 
of the type 

e{no) e{no) 

Hp{z,no) = co,+(3{no)z-P-+^ J] (z - Mno))"'^"''\ ^ Sk{no) = p, 

k=l k=l 

i^kino) + vw [no] for fc fc', fc, fc' = 1, . . . , £(no), z £ C \ {0}, (4.32) 
where we introduced the coefficient /3{no). We define 

Vk{no) = {vkino), ~{2/co^+)vk(no)^- Gp{vk[no),no)), fc = 1, . . . ,£(no). (4.33) 
An explicit computation of /3(no) then yields 

,(no) (d^''=(-)-i(C-iy(P))/dC^'=("'''-i)|p^(^,^)^^^(„„) 



a+(no)/3(no) = -\^ 



2 (Pk{no) - 1) 
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9(no) 

k' = l, k'=ik 

q{no) 2p+l 



+ K^Sf ^ ^ - > >™ )• (4.34) 



fe=l m=0 ^ 



The result (|4.34p is obtained by inserting the expressions (|4.18p . (|4.19p . and (|4.32p 
for • , no), Gp{- ,no), and Hp{- ,no) into (|4.3ip and collecting all terms of order 

(VI) Introduce 

^(no) /2p+l \ -1/2 

P+{no) = /3(no) n '^^("0)^''^"°^ ( n ^"0 • (4.35) 

^ m=0 ' 

(VII) Using Gp(z, no), Hp{z,no), Fp{z, hq), /3(no), Q;+(no), and /?+(no), we next 
construct the no ± 1 terms from the following equations: 

F- = —{a^Hp - 2aGp + Fp), (4.36) 
— Z7 — — - 



z 



respectively, 



H- = -{13' Fp - 2[3Gp + Hp), (4.37) 

^ - - - 

GI ^-{{l + a(3)Gp-aHp- (3Fp), (4.38) 
- 7 _ _ _ 



F+ = -l-{{a+)^Hp + 2a+zGp + z^Fp), (4.39) 

H+ = -!-{{(]+ z)'Fp + 2p+zGp + Hp), (4.40) 

— — — — 

G+ = ^((1 + a+P+)zGp + a+Hp + P+z^Fp). (4.41) 

- 27+ £1 £1 



Moreover, 



(Gp )^ - ^^^-i/p = Rp, (G+)2 - - (4.42) 

Inserting (gUl), gH), and (l432ll in (|436l)-(|438l) one verifies 

F-(z,no) - -co,+a(no)^^+"' + 0(zP+-2), (4.43) 

— Z^OO 

H-{z,no) - 0(zP+), (4.44) 
F-(z,no) ^0{z-P-), (4.45) 

— 

i?„~(z,no) = -co,-l3{no)z-P-+^ +0{z~P-+^), (4.46) 

G^(z,no) = ico.-z-^"- +--- + ico,+0''+. (4.47) 
The last equation implies 

Gp(z,no-l)'-i?p(^) = 0{z'P+-^), (4.48) 

— — ^^00 

Gp{z, no - 1)2 - Rp{z) 0{z-^P-+^), (4.49) 

so we may write 

Gp(z,no-l)2-i?p(z) = F,,p^_i(z,no-l)^p__i,.(2,no-l), z e C\{0}, (4.50) 
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for some s G r G where 

^s,p+-l("-o - 1) = C-sZ^" H Co,+Q:(no)2:^+"\ 

Hp_-i,rin,o - 1) = -co,-/3(no)z"^"+^ H h c^z''. 

The right-hand side of (|4.50p cannot vanish identicaUy (since otherwise Rp{z) 
Gp{z, no — 1)^ would yield double zeros of Rp{z)), and hence, 



no - 1) 



Gp{z,no - 1) + {co^+/2)z P-y _ i?p__i,r(z, no - 1) 



Fs,p_^_-i{z,nn - 1) Gp(z,no - 1) - (co,+/2)z-P-y ' 

P=iz,y)elCp. (4.51) 

Next, we decompose 

P+-1+S 

^'s,p+-i(^,?^o - 1) -co,+a(no)2;"'' J]^ (2; - /ij(no - 1)), (4-52) 

p_ — 1+r 

7?p__i,,(z,no-l) = C^^"''-+' n (4.53) 

where C G C \ {0} and {iij{no - 1)}^;7^+' C C, {i/j(no - l)}^!^^^'' C C. The 
divisor of 0( • , no ^ 1) is then given by 

((^(•,no-l)) = 2?Po,_&(„o-i) - 2?p<^_A("o-i)' (4-54) 

where 

p.{no - 1) = {^'o.+, ■ • . ,Po,+ , Ai(»^o - 1), ■ ■ ■ ,Ap+-i+s('^o - 1)}, (4.55) 

~ " v ' 

P— — s times 

E.{no - 1) = {i>i(no - 1), . . . , i>p__i+r(no - l),Poo+, • ■ ■ ,-Poo+}- (4.56) 

—r times 

In particular, 

^/i(no-i) is an admissible divisor if and only if s = p-, (4-57) 

I?£,(„P_i) is an admissible divisor if and only if r = p^. (4.58) 

(VIII) Assuming that ()4.48p . ()4.49p are precisely of order 2;±(2p±-i)^ that is, as- 
suming s — p^ and r = in ()4.5H) , we rewrite (|4.51|) as 

Gp(z,no-l)2-i?p(z) = i^p(z,no-l)i?p(z,no-l), zgC\{0}. (4.59) 

By construction, i^p( • , no — 1) and Hp{ • , no — 1) are then of the type 

g(rio-l) 

Fp(z,no-l) = -co.+a(no)z-f- \[ (z - /i,(no - l))^"=("°-i), 

ij(no-l) 

Pk{m-l)=P, (4.60) 

k=l 

A*fe('^o - 1) ^ ^J■k'{na - 1) for k ^ k' , k,k' = 1,.. . , g(no - 1), z G C \ {0}, 

e{no-l) 

Hp{z,no-l) = co.+l3{no~l)z-P-+^ [] _ _ l))-.("o-i)^ 

fc=i 



22 F. GESZTESY, H. HOLDEN, J. MICHOR, AND G. TESCHL 

t(no-l) 

^ Sfc(no-l)=p, (4.61) 
fe=i 

Vk{no-l)^Vk'{no-l) for k ^ k' , A:, fc' = 1, . . . , ^(no - 1), zeC\{0}, 

where we introduced the coefficient /3(no — !)• We define 

A/c("o - 1) = (Mfc("-o - 1), (2/co,+)Affe("-o - l)^- Gpinkino - l),7io - 1)), 

fc=l,...,g(n„-l), ^^^^^^ 
:^fc(no - 1) = (i^fe(no - 1), -(2/co,+)i^fe(no - 1)^ Gp{i^k{no - l),7io - 1)), 

fc = l,...,^(no-l). 

(IX) At this point one can iterate the procedure step by step to construct Fp{- ,n), 
Gp{ ■ , n), Hp{ ■ , n), a(n), /3(n), Hj{n), Vj{n), etc., for n € (— oo, no] fl Z, subject to 
the following assumption (cf. (|4.57p . (|4.58p ) at each step: 

admissible divisor (and hence a{n — 1) 7^ 0) (4.63) 
for all n e (—00, rto] n Z, 
T^i>{n-i) is an admissible divisor (and hence /3(n — 1) 7^ 0) (4.64) 
for all n G (—00, tiq] H Z. 

The formalism is symmetric with respect to tiq and can equally well be developed 
for n G (—00, no] n Z subject to the analogous assumption 

T^ii(n+i) is an admissible divisor (and hence a{n + 2) 7^ 0) (4.65) 

for all n e [no, 00) n Z, 

■^^(n+i) admissible divisor (and hence /3(n + 2) 7^ 0) (4.66) 

for all n G [no, 00) n Z. 

(X) Choosing the initial data such that 

%„o) G -^0, (4.67) 

where AIq C Sym^(/Cp) is the set of admissible initial divisors introduced in Lemma 
14.31 then guarantees that assumptions (|4.63p - (l4.66p are satisfied for all n G Z. 

Equations (|4.36p - (|4.4ip (for arbitrary n G Z) are equivalent to s-ALp(a, (3) = 0. 

At this stage we have verified the basic hypotheses of Section [3] (i.e., (|3.2p and 
the assumption that a,/? satisfy the pth stationary AL system (j2.33p ) and hence 
all results of Section [3] apply. 

In summary, we proved the following result: 

Theorem 4.4. Let n G Z, suppose the set {Em}m=o satisfies E„i 7^ for 
m 7^ m', m,m' — 0, . . . ,2p + 1, and introduce the function Rp and the hyperelliptic 
curve fCp as in (|3.ip . Choose a{no) G C \ {0}, co.+ G C \ {0}, and a nonspecial 
divisor Vj^i^l^na) ^ -^o, where Aio C Sym^(A^p) is the set of admissible initial divisors 
introduced in Lemma [4.31 Then the stationary {complex) Ablowitz-Ladik algorithm 
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as outlined in steps (I)-{X) produces solutions a, (3 of the pth stationary Ablowitz- 
Ladik system, 

^ V/^(V.+ +5p-.-) + - V-i.- / (4.68) 



P={P-,P+) e 



n^"0 ^(n,no)a(no), (4.69) 



satisfying (13. 2p and 

•2p+l s (n-no)/2 



m— 



9(n) (rfP;=(")-l {C^y{P)) /dCP''(")-l) 



J'=(C,';)=/ifc(") 



2 iPkin) - 1)! 

g(n) 

X n (Mfc(")-A^fe'(ri))-'"'^") 
fc'=i, fe'^^fc 

- 2p+l . 1/2 ?(") 



, .2p+l .1/2 g(n) 9(n) 2p+l . \ 

^ \ m=0 ^ fc=l fc=l m=0 ^ / 
2p+l X -(n+l-no)/2 

\{E„A A{n + l,noy^a{noy\ (4.70) 



where 

y4^(n, no) = < 1, ?^ = ^^o, (4-71) 

WClYlt,^ l^k{n'y^^-'\ n<no-l. 
Moreover, Lemmas I3.2H3.5I apply. 

Finally, we briefly illustrate some aspects of this analysis in the special case 
p — (1, 1) (i.e., the case where p.ip represents an elliptic Riemann surface) in more 
detail. 

Example 4.5. The case p = (1, 1). 
In this case one has 

F(iA){z,n) = -co,+a(n + l)2"^(z - ^i(n)), 

G(M)(^>'^) = ^(^-^^(;^)^(i,i)(^,«) + ^(i,i)(Ai("))'/', (4.72) 

and hence a straightforward calculation shows that 

G(i4)(z,n)^ - i?(i^i)(z) = -Co^^Q:(n + l)/3(n)z~^(z ~ ^i(n))(z - i'i(n)) 

y(Ai(n)) , ^1/2 



c. 



{z - fii (n)) — \ — + ^1 (n) — z 



2z \\ tJ-i{n) Mi('^) 2 



(4.73) 
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E ( 1 y{Mn)) , MiH , 1 



MiW V E^l'^ A<i(n) £;i/2 ^i(n) 2 
w/iere 

^^=E^™' ^=ni?™. (4.74) 

m— m— 

Solving for i'i{n) one then obtains 

Ml(?l) ^ I ^^(n) I 1 E- ^ ' 

Mi(n) £1/2 ^ ^i(n) 2 

Thus, z/i(no) could be or oo even i//ii(7io) ^ 0, cxd. 

5. Properties of Algebro-Geometric Solutions 
OF THE Time-Dependent Ablowitz-Ladik Hierarchy 

In this section we present a quick review of properties of algebro-geometric so- 
lutions of the time-dependent Ablowitz-Ladik hierarchy. Again we omit all proofs 
and refer to [29j, [31], and [32] for details. 

For most of this section we assume the following hypothesis. 

Hypothesis 5.1. (i) Suppose that a, P satisfy 

a{- ,t),l3{- ,t) eC^, te M, a{n, •), f3{n, •) e C^(R), n G Z, 

a{n, t)P{n, t) ^ {0, 1}, (n, t) £ Z x M. ^^'"^^ 

(m) Assume that the hyperelliptic curve ICp satisfies p.l|) and ([33]). 

In order to briefly analyze algebro-geometric solutions of the time-dependent 
Ablowitz-Ladik hierarchy we proceed as follows. Given p G Ng, consider a complex- 
valued solution q;(°) , /S^^^ of the pth stationary Ablowitz-Ladik system s-ALp(a, b) = 
0, associated with ]Cp and a given set of summation constants {c£^±}^=i_...^p^ C C. 
Next, let r G Ng; we intend to consider solutions a = a{tr),P = P{tr) of the rth 
AL flow AL^(a,/3) = with a(io,r) = a^"\/3{to,r) = /S^"^ for some ig,^ G M. To 
emphasize that the summation constants in the definitions of the stationary and 
the time-dependent Ablowitz-Ladik equations are independent of each other, we 
indicate this by adding a tilde on all the time-dependent quantities. Hence we 
shall employ the notation Vr, Fr, Gr, Hr, K^, fs.±, gs,±, ^s,±i Cs,±, in order to 
distinguish them from Vp, Fp, Gp, Hp, Kp, fi^±, gi,±, hi^±, ci^±, in the following. In 
addition, we will follow a more elaborate notation inspired by Hirota's r-function 
approach and indicate the individual rth Ablowitz-Ladik flow by a separate time 
variable tr G M. More precisely, we will review properties of solutions a, (3 of the 
time-dependent algebro-geometric initial value problem 

( n\ (-^'^t^~'^i9r+.+ +gr^,-)+ fr+-l.+ - fr_-iA „ 

y-i/3t,+^(5r+,+ +5r-,-)-^r_-i.- + ^.^_i,+y (5.2) 
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for some to,r G M, where a = a(n,tr), /3 = P{n,tr) satisfy (|5.ip and a fixed curve 
ICp is associated with the stationary solutions a^'^\ P'-^'> in (|5.3[) . Here, 

£=(P-,P+)eN2\{(0,0)}, 1= (r_,r+) eN^, p = p_+p+-f. (5.4) 

In terms of the zero-curvature formulation this amounts to solving 

Ut^{z, tr) + U{Z, U)Vr{z, tr) - V+{z, tr)U{z, U) = 0, (5.5) 
U{z,to,r)Vp{z,t^.r) - V+(z,tQ,r)U{z,ta,r} = 0. (5.6) 

One can show (cf. Lemma that the stationary Ablowitz-Ladik system (|5.6p is 
actually satisfied for all times tr_ € R. Thus, wc impose 

Ut^ + UVr^-V+U = 0, (5.7) 

UVp-V+U ^0, (5.8) 

instead of (15. 5|) and ()5.6p . For further reference, we recall the relevant quantities 

here (cf. dMOD, (E^D, UM-^M, UM)- 



U{z)- 
Vp{z)^i 



z a 
zfi 1/' 



^Gp(z) ~F^{zf 
H-{z) -G-{z) 



Vr{z) = i 



,(G-{z) ~-F-{zy 



(5.9) 



^H-{z) -K-{z) 



and 



Fp{z 
Gp{z 
Hp{z 

Fr_{z 
Gr{z 

Hr{z 
Kr{z 



^fp--£-z ^+ fp^-i-e,+z'^ ^ -co,+a+z p Yl{z - ^J.J), 

1=1 1=0 j = l 

p- p+ 

Ygp^-i,+z\ 



Ygp--L-z~ 
1=1 



l=Q 



1=0 



1=1 j=i 

r_ r+ — 1 

s=l s=0 
r_ r-i- 

Ygr_-s-z~'' + Ygr+-s,+z'', 



(5.10) 



s=l 



s=0 



s=0 



E 



gr_-s~Z " + 



»■+—«>+ 



= Gr(2^) + .9r_,- — g 



for fixed £ G \ {(0,0 )}, r E N^. Here fe^±, fs,±, ge,±, gs,±, he,±, and hs^± are 
defined as in (I2.5p - (l2.12p with appropriate sets of summation constants Ci^±, I G No, 
and Cfe±, k E Nq. Exphcitly, ((^ and (jO)) are equivalent to (cf. 



(5.11) 
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Pu — 






= 


z{G:; - Gr) + zl3Fr + aH^ , 


(5.13) 


= 


ZI3F.;' + aHr^ + K:; - Kr^, 


(5.14) 


= 


z(G- - Gp) + zpFp + aH- , 


(5.15) 


= 


zPFp + aHp -Gp + Gp, 


(5.16) 


= 


-Fp + zFp +a{Gp + Gp), 


(5.17) 


= 


zP{Gp + Gp) - zHp + Hp, 


(5.18) 



respectively. In particular, (j2.39p holds in the present tr-dependent setting, that is, 

Gl~FpHp^Rp. (5.19) 
As in the stationary context (|3.8p . (13. 9[) we introduce 

p,j{n,tr) = (/ij(n,ir), (2/co,+)Aij(n,tr)^"Gp(^j(n,t^),n, ir)) e /Cp, 

j = 1, . . . ,p, {n,tr) e Z X M, ^ ' ^ 

and 

Vj(n,tr) = (^■j(n,ir), -(2/co,+)^'j(n,tr)''"Gp(^'j(n,<r),"-,ir)) G ^p, 

] = l,...,p, (n,tj e Z X M, ^^'^^^ 

and note that the regularity assumptions (|5.ip on a,/3 imply continuity of ^.j and 
i^fe with respect to G K (away from collisions of these zeros, and Vk are of 
course C°°). 

In analogy to (|3.12p . (|3.13p . one defines the following meromorphic function 
4>{ ■ , n,tr) on ICp, 

^, (co,+/2)z-P-y + Gp(z,n,y 
</'(P, i.) = F(znt) ^^-^^^ 

-Hp_iz,n,U) ^^^^^^ 



(co,+/2)z P-y - Gp{z,n,tr)' 
P = (z,y) e /Cp, (n,tr) e Z X M, 
with divisor (0( ■ ,n,tr)) of (/)( • , 7i, i,-) given by 

(0( • , tr)) = 'Dpo,-P(n,U) - ^Poo_ A(",tr:)- (^-2^) 

The time-dependent Baker-Akhiezer vector is then defined in terms of (j) by 

*,p,„,„.,.,M^(*;j^:::::-::j). (5-^^) 

V'i(-P, n,no,ir, io,r) = exp (^i^ ds(Gr(2;, no, s) - Fr(2;, no, s)0(P, no, s))^ (5.26) 

nI!'=„o + l + y ^rj) , n > 710 + 1, 

1, n = no, 

.n"L„+i + a(n',y0-(P,n',t,))"\ n < no - 1, 

tp2{P,n,no,tr,ta,r) = cxp (^i ^ ds(Gr(2:, no, s) - f;(z, no, s)0(P, no, s))^ (5.27) 



X < 
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f n^'=„o+i ^l)-^" (P, n', y -1 + 1) , n > no + 1, 

X (/)(P, no,tr) < 1, 71 = 710, 

in"U+i (z/3(7i',i,)0-(F,n',M-i + iy\ 11 < no - 1, 
P = (z, y)elCp\ {Poo+ , Poo_ , Po.+, Po,-}, {n, U) E Z x K. 

One observes that 

7/'l(P, 71, 710, tr, ir) = IpliP, "o, TlQ, U, ir)lpi{P, Tl, rio^U, U) , 



P = {z,y) e ICp\ {Poo+,Poo_,-Po,+ ,-Po,-}, {n,nQ,tr,ir) ^ 



X 



(5.28) 



The following lemma records basic properties of and ^' in analogy to the 
stationary case discussed in Lemma [321 



Lemma 5.2 ([31J. Assume Hypothesis \5.1\ and suppose that ()5.7|) . (|5.8p hold. In 

addition, let P = {z,y) G ICp\ {Poo+,^oo_}, (n, tio, ^r, iQ,r) G x M^. T/ieTi (/) 
satisfies 

a4>{P)4>- (P) - (P) + z(/)(P) = z/3, (5.29) 

(j)tAP) = iF,4\P) - i(G.(z) + ^,(z))(/.(P) + iH,(z), (5.30) 

mHP*)^-^, (5.31) 
Gp{z) 

</.(P) + </.(P*) = 2-^, (5.32) 

</.(P)-</.(P*) = co.+z-^'-|^. (5.33) 

Moreover, assuming P — [z, y) e /Cp \ {Poo+, Poo_ , -fo,+ , -Fb,-}; ^''■e'^ ^ satisfies 

■02(^,71, 7io,tr,io,r) = (l){P,n,tr)il^i{P,n,no,tr,to,r), (5.34) 
t/(z)*-(P) = *(P), (5.35) 
Fp(z)*-(P) - -(i/2)co,+z-f-y*-(P), (5.36) 

^tAP) = V+{z)-^iP), (5.37) 

Fp{z,n,tr) 

ipi{P,n,no,tr,to^r)i'i{P*,n,no,tr,tQ,r) = z"' ""-7^ r-r(71,71o,tr), (5.38) 

- - - - Pp(z,7lo,io,r) 

Hp{z,n,tr) 

ip2{P,n,nQ,tr,tf)^r)'4>2{P* ,n,no,tr,tQ,r) = z^ ""tt-t : — ^^{n,na,tr), (5.39) 

- - - - Pp(z,7lo,io,r) 

ipi (P, 71, TiQ, tr, io,r)V'2 {P* ,n,no, tr, to.r) + Tpl (P* , 71, no, t^, to,£)V'2 [P, n, no, tr , tox) 
Gr) [z, n, tr) 

= 2z"-""— P ^r(n,no,t.), (5.40) 

Pp(z,no,io,r) 

■01 {P, n, no, tr, to,r)^2 {P* ,n,no, tr, to,r) - ipi [P* , n, no, tr, to,r)ip2 [P, n, no, tr, to,r) 

= -co,+z^-^°-P- --^_-r(n, no, tr), (5.41) 
rp[z, no, iQ r) 
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where 

{nl^'=«o+i'T("',ir) n>no + l, 
1 n = no, (5.42) 

n:L„+i7«y-^ n<no-l. 
In addition, as long as the zeros ijLj{nQ,s) of {■)p- Fp{- ,no, s) are all simple and 
distinct from zero for s G Ip,, 1-^ ^ M an open interval, ^(-,71, tiq, t^-, to,r) is vfiero- 
morphic on JCp \ {Poa+ , ^'oo- , -Po,+ , -Po,- } for {n, tr, to,r) xlf^. 

The trace formulas recorded in Lemma |3 . 3 1 extend to the present time-dependent 
context without any change as G M can be viewed as a fixed parameter. Further 
details are thus omitted. 

For completeness we next mention the Dubrovin-type equations for the time 
variation of the zeros /ij of {■)P^Fp and lyj of {■)P-~^Hp governed by the AL^ 
flow. 

Lemma 5.3 f[31j). Assume Hypothesis 15.11 and suppose that (j5.7p . (|5.8p hold on 
Z X Xp with Xp C W an open interval. In addition, assume that the zeros fij, j — 
1, . . . ,p, of {■ Fp{ ■ ) remain distinct and nonzero on Z x X^. Then {fj-j}j=i,...,p, 
defined in (|5.20p . satisfies the following first- order system of differential equations 
on "Z X Ip, 

^ p 
^iJ,u = -iFr{^^ij)y{ilj){a+)-^Y\_^l^3 ~ l^k)~^^ J = l,---,?5, (5-43) 

fc=i 

with 

(i,{n,-)&C^{lM^ j = nGZ. (5.44) 

For the zeros Vj , j = 1, . . . ,p, of {■ Y^^^Hp{ ■ ), identical statements hold with 
and Jp replaced by Vj and T^, etc. {with Xiy C_ R an open interval). In particular, 
{£'j}j=i,...,p, defined in ()5.2ip . satisfies the first-order system on Z x X^, 



iHAyj)y(yj){^Vj) ^^{v^-Vk) \ j^l,...,p, (5.45) 



fc=i 



with 



i>j(n,-) G C°°(J^,/Cp), j = nGZ. (5.46) 

When attempting to solve the Dubrovin-type systems (|5.43|) and (|5.45[) . they 
must be augmented with appropriate divisors X>fi(noM r) ^ Sym^ /Cp, to,r G 2^, and 
T^o{no.ta r-) £ Sym^ ICp, to,r € Xi,, as initial conditions. 

For the ^r-dependence of Fp, Gp, and Hp one obtains the following result. 

Lemma 5.4 ([31 ). Assume Hvvothesis \5.1\ and suppose that (|5.7p . (j5.8p hold. In 
addition, let {z,n,tr) G C x Z x M. Then, 

Fp^t^ = -2iGpFr + i (G^ + K^r) Fp, (5.47) 
Gp^t^ = iFpHr - iHpFr, (5.48) 
Hp^t^ = 2iGpHr~i{Gr + Kr)Hp. (5.49) 
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In particular, (I5.47P - (|5.49|) are equivalent to 

Vp^t^ = [V;, Vp] . (5.50) 

It will be shown in Section [S] that Lemma 15.41 yields a first-order system of 
differential equations for fi^±, gi.±, and hi^±, that serves as a pertinent substitute 
for the Dubrovin equations ()5.43p even (in fact, especially) when some of the /ij 
coincide. 

Lemma [3. 5 1 on nonspecial divisors and the linearization property of the Abel map 
extend to the present time-dependent setting. For this fact we need to introduce 
a particular differential of the second kind, ftr , defined as follows. Let ujp 

— J oo _|_ 7 y 

(2) 

and _^ q be the normalized differentials of the second kind with a unique pole 
at Poo± and Po,±7 respectively, and principal parts 

^pI^,g^=o (r'"' + 0(l))rfC, P^Poo^, C^l/z, qeNo, (5.51) 



OJ 



(2) 

Po,±,q 



(C"'"'^ + 0(l))dC, P ^ ^o,±, C = ^, 9 e No, (5.52) 



with vanishing a-periods. 



Moreover, we define 

o(2) - ^ fsr (, ,(2) , ,(2) 



Pn,-,s-i; 
s=l 

'+ \ 



(5.54) 

s=l ' 

where q.± are the summation constants in F^. The corresponding vector of h- 
periods of £7^ /{^-nii) is then denoted by 

= (C/(^,\...,[/(^)), m^^l^^, J = (5.55) 



The time-dependent analog of Lemma [3T5l then reads as follows. 

Lemma 5.5 f[31|. [32]). Assume Hvvothesis \5.l\ and suppose that (j5.7[) . (|5.8p /loM. 
Moreover, let {n,tr) S Z x R. Denote by Vfi, fi — {fii, . . . , fip} and Vc^, p_ — 
{vi, . . . ,i'p}, the pole and zero divisors of degree p, respectively, associated with a, 
(3, and (f> defined according to ()5.20l) and (|5.21|) . that is, 

flj{n,tr) = {flj{n,tr),{2/cQ,+ )flj{n,trf-Gp{^lj{n,tr),n,tr)), j = l,...,p, 

(5.56) 

i>j{n,tr) = {iyj{n,tr),-{2/co^+)vj{n,try^Gp{vj{Ti,tr),n,tr)), j = 1, 

(5.57) 

Then 'Dj^^(^n.tr) o.'^d 'Du(n,tr) '^'"^ nonspecial for all {n,tr) G Z x R. Moreover, the 
Abel map linearizes the auxiliary divisors 'Dp(n.tr) '^'^^ T^i>{n,tr) *^ sense that 

~ (2) 

^Qo = ^Qo (%no,to,.)) + no)Apo^_ (^oo+ ) + {tr- to,r)IZr > (5-58) 
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^Qo = aQo(%("o,to..)) + (" - no)Ap^^ _ (Poo+) + {tr - to,r)Ur , (5-59) 

~(2) 

where Qq G /Cp is a given base point and [/,. is the vector of b-periods introduced 
in ([^35)1 . 

Again we refer to [35] (and the references cited therein) for theta function rep- 
resentations of a and f3. These theta function representations also show that 
'y(n,tr) ^ {0,1} for all {n,tr) £ Z x R, that is, the last condition in (|5.ip is 
satisfied for the time-dependent algcbro-geomctric AL solutions discussed in this 
section provided the associated Dirichlet and Neumann divisors are admissible. 

6. The Algebro-Geometric Ablowitz-Ladik Hierarchy 
Initial Value Problem 

In this section we consider the algebro-gcometric Ablowitz-Ladik hierarchy ini- 
tial value problem ()5.2|) . ()5.3p with complex- valued initial data. For a generic set 
of initial data we will prove unique solvability of the initial value problem globally 
in time. 

As mentioned in the introduction, the strategy of basing the solution of the 
algebro-geometric initial value problem on the Dubrovin-type equations (j5.43p . 
()5.45p . and the trace formulas for a and (3, meets with serious difficulties as the 
Dirichlet eigenvalues fij may collide on ICp and hence the denominator of (|5.43|) can 
blow up. Hence, we will develop an alternative strategy based on the use of ele- 
mentary symmetric functions of the variables {fJ'j}j=i,...,p in this section, which can 
accommodate collisions of fij. In short, our strategy will consist of the following: 

(z) Replace the first-order autonomous Dubrovin-type system (|5.43p of differen- 
tial equations in tr for the Dirichlet eigenvalues fij{n,tr), j = 1, . . . ,p, augmented 
by appropriate initial conditions, by the first-order autonomous system (|6.57p . 
(|6.58p for the coefficients fi,±, /i£,±, i = 1, . . . ,p± — 1, and gi,±, £ = 1, . . . ,p±, 
with respect to tr- Solve this first-order autonomous system in some time interval 
{to.r — 7o,io,r + Tq) under appropriate initial conditions at (no,io.r) derived from 
an initial (nonspecial) Dirichlet divisor 'Djx^no.ta r)- 

{a) Use the stationary algorithm derived in Section 2] to extend the solution of 
step {i) from {tiq} x (to.r — To, io.r + To) to Z x (to.r — Tq, to.r + Tq) (cf. Lemma 

(Hi) Prove consistency of this approach, that is, show that the discrete algo- 
rithm of Section [4] is compatible with the time-dependent Lax and zero-curvature 
equations in the sense that first solving the autonomous system (I6.57p . (I6.58P and 
then applying the discrete algorithm, or first applying the discrete algorithm and 
then solving the autonomous system (|6.57p . (|6.58p yields the same result whenever 
the same endpoint {n,tr) is reached (cf. Lemma [6731 and Theorem 16. 4p . 

(iv) Prove that there is a dense set of initial conditions of full measure for which 
this strategy yields global solutions of the algebro-geometric Ablowitz-Ladik hier- 
archy initial value problem. 

To set up this formalism we need some preparations. From the outset we make 
the following assumption. 

Hypothesis 6.1. Suppose that 

a,(3 eC^ and a{n)l3{n) ^ {0, 1} for all n e Z, (6.1) 
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and assume that a, (3 satisfy the pth stationary Ablowitz-Ladik equation (|2.33p . 
In addition, suppose that the affine part of the hyperelliptic curve ICp in (13. ip is 
nonsingular. 

We introduce a deformation (time) parameter t^. £ M in a = a{tr) and (3 — (3{tr) 
and hence obtain t^-dependent quantities fi = fi{tr), gg = gi{tr), Fp{z) = Fp{z, t^, 
Gp{z) — Gp{z,tr), etc. At a fixed initial time to,r S we require that 

(a,/?)|t^.,„.^ = /3(")), (6.2) 

where a^°^ — a{ - , io.r)i Z?''*''' = /3( • , ^o.r) satisfy the pth stationary Ablowitz-Ladik 
system (|2.33p . As discussed in Section |31 in order to guarantee that the stationary 
solutions (|6.2p can be constructed for all n E Z one starts from a particular divisor 

%„o.*o,J e Mo, (6.3) 

where fl(no,to,r) is of the form 

A("-o,^o,r) (6.4) 

= {ft-lino, to,r), . . . , fil{no,to,r), ■ ■ ■ , Ag(no,to (^^0 , ^0,r) , ■ • • , Ag(no.to r ) ('^0 , ^0,r ) } ■ 
. ' ^ — ' 

Pl(nO:to.r) times Pg(no ,to.r ) ("0 .*0,r ) times 

Moreover, as in Section 0] we prescribe the data 

a{no, io,r) e C \ {0} and co,+ G C \ {0}, (6.5) 

and of course the hyperelliptic curve ICp with nonsingular afRne part (cf. (|4.10p V 
In addition, we introduce 

. 2p+l . 1/2 g(«o,*o,r) 

a+{no,to,r) = a{no,to,r)[l[E„A [J A^fc(^o, M"'"''""'*"'^^ (6-6) 

^ ?ri=0 ^ fc=l 

Fp(2:,no, io,r) = ^ /p_-f -(»T-o, io,r)2"^ + ^ fp+-i-e,+ (na,to,r)z^ 
i=\ 1=0 

= -co^+a+{no,to,r)z-P- J] " A^^K, (6.7) 

k=l 

Gp{z,no,to,r) = ^( — r~7 — T +r~l — AFp{z,no,to.r) 



z 



-Fp(z,no,io,r) 



2a+{no,to,r) 

q(no,to,r) Pk{no,ta.r)-l (d^ (^^^y (P)) / d^^ 



t[ U £\{pk{no,to,r)~i~l)l 



X 



«("0,*0,r) ^ 

X n (C-A'fc'(^io,io,r))-^'='("°^*"-' 

k' = l,k'^k / 

in analogy to (|4.16p . 



C=/ifc(n0,t0,r;) 
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Our aim is to find an autonomous first-order system of ordinary differential 
equations witli respect to tr for fi^±, gi,±, and hi_± rather than for /ij. We divide 
the differential equation 

^ -2iGpF,_ + i{Gr + kr)Fp (6.9) 



by Co, +2: P-y and rewrite it using Theorem lC.il as 

e=o ^ s=i s=o e=o 

+ z ( 2 ^ + 2 ^ 5r + -«.+C^' - ffr-,- + K,+ ) E A^+C'+'. (6-10) 

^ s=0 s=l ^ £=0 

P^P^^, C = l/z. 

The coefficients of C~*, s = 0, . . . , r+ — 1, cancel since 

E h-k,+h,+ = E 9£-k,+fk,+, t e No. (6.11) 

fc=0 fc=0 

In ([6TT|) we used ([2T9l) . 

Comparing coefficients in (I6.10p then yields Q 

'■+-1 

f£.+,u = ife,+ igr+.+ - 5r_ -) + 2j E {9k.+fr++e-k,+ - /fc,+5r++«-fc,+) (6.13) 

e e 

-2i ^ gk,+fr_~i-£+k,~ + 2i ^ fk^+gr_-£+k,-, ^ G No. 

fc=(£+l-r_)VO fc=(£+2-r_)V0 

By (|6.1ip . the last sum in (|6.13p can be rewritten as 

r+-l 

E (ffi: + /r++«-j,+ - f],+gr++£-j,+) 

/ r++£ r++£\ 
= ( E ^ E ) {93, + Ir++£-0,+ - /j,+ffr++£-j,+) 

= ^ E ^ f].^+gr++l-], + ) 

e. 

= E + - /j,+5r++£-j, + )- (6.14) 

One performs a similar computation for fe.~,tr using Theorem lC.il at P Pq,+ - 
In summary, since fk,± — '^£—QCk—i,±f£,±j l|6.13|l and (16.1411 yield the following 



*m V n = max{m, n}. 
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autonomous first-order system (for fixed n ~ uq) 

fe,±,tr,^ ^e,±{fj-,fj,+ ,9j-,9j,+), £ = 0, . . . ,p± - I, (6.15) 
witli initial conditions 

fe,±{m,to^r), ^ = 0, . . . ,p± - 1, (6.16) 
5f,±(»^o,to,£), £ = 0,...,p±, 
wfiere J-e.±, i — 0, . . . ,p± — 1, are polynomials in 2p + 3 variables, 

I 

+ $Z (/fc,±(5r^-«+fc,=F ~ 9r±+l-k,±) + 3fc,± (/r±+f-fe,± " /r:p - l-^+fc,=F )) 
fe=0 

I fO, < fc < - 1, 

[ - Ei=0~''^ fj,±9r^-k+j,T^ fc > rzp - 1. 

Explicitly, one obtains (for simplicity, r± > 1) 

^n,± = ifo,±i9r^:T - 9r±.±) + 2i50,±(/r±.± " (6.18) 

^i,± = 2i/o,±(gr:p-i,T - 5r±+i,±) + i/i,±(ffr:p,T - 9r±.±) 

+ 2i50,±(/r±+l,± " fr^-2,T) + 2jgi,±(/r±,± " /r:p-l,T)' ^tC. 

By ()6.6p ^ (|6.8p . the initial conditions (|6.16p are uniquely determined by the initial 
divisor I'/i(„g,fo ^) in (|6.3p and by the data in (|6.5p . 
Similarly, one transforms 

Gp,*^ = iFpi/^ - iHpFr^, (6.19) 

Hp = 2iGpHr - i{Gr + Kr)Hp (6.20) 

into (for fixed n = no) 

5o,±,t, = 0, 

r± — 1 

9£,±,t^ — i ^ {hk.,±fr±-l+t-k,± - /fe,±/lr±-l+£-fc,±) 
fe=0 

fc=(f-r^)VO 
£-1 

— — fk,±hr±-l+l-k,±) 

k=0 

t-l 

+ i ^ (/fe,±^r:p-£+fe,T - ^fe,±/rT-''+fc,T)' ^ £ N, (6.21) 
/c=(£-r^)VO 

r± — 1 
k=0 



y n = max{m, n}. 
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I 

= iht^±[gr^,^ — 3r±,±) + (^fe,±3r±+£-fe,± — 9k,±K±+e-k,±) 

k=0 

I I 

fe=(£-r^ + l)VO fc=(£-r^)VO 

^ e No. (6.22) 

Summing over I in (|6.2ip , (|6.22p then yields the foUowing first-order system 

g(,±,U=Ge,±{fk-,fk,+,hk-,hk,+ ), i = 0,...,p±, (6.23) 
he,±,tr_ = 'H£^±{gk,-,gk,+ ,hk-,hk,+), £ = 0, . . . ,p± - 1, (6.24) 
with initial conditions 

/«,±('^0,i0,r), ^ = 0, . . . ,p± - 1, 

9e,±{no,to^r), £ = 0,...,p±, (6.25) 
hi,±{no,to^r), ^ = 0, . . . ,p± - 1, 
where ^f,±, 'He,±, are polynomials in 2p + 2, 2p + 3 variables 

£-1 

Gl,± — «^ (/fc,±(/lr^-f+fc,ip — hr±-l+l-k,±) + hk^±{fr±-l+e-k,± — fr^-l+k,^)) 



k=0 
-1 



^ 0, < fc < r 

-l2_^Ci-i^k,±X< k-r^-l,f 7 f ? \ l.^ 



(6.26) 

+ 2i^ (gfc,±(^r:p-l-£+fc,T ~ ^r±+£-fc,±) + ^fc,±(5r±+£-fc,± " 5r^ +fc,T )) 
A:=0 

'O, 0<fc<rzp, 
+ 2i ^ ci-k,± X - Ej=o'' 5j\±^r^-i-fc+j,T (6-27) 



fc=0 



Explicitly (assuming r± > 2), 
Go.± - 0, 

0i.± = i.fa,±{hr^-i,T - hr±,±) + i/io,±(/r±,± - fr^-l,^)^ (6-28) 

Q2,± — */o,±(/lr:p-2,=F ~ hr± + l,±) + */l,± (/ir:^ -1,=F ~ hr±,±) 

+ i/jO,±(/r± + l,± - fr^-2,T) + *^l,±(/r±,± " etC, 
Ho,± 2i5ro,±(/lr:p-l,T - ^r±,±) + i/l0,±(5r±,± " ffr:p,T)> (6.29) 
Wl^i = 2igo,±(/lr:p-2,T ~ hr^ + l^±) + 2igi^±(/lr^_i^zp ~ ^r±,±) 

+ 2i/lo,±(gr± + l,± - gr:p-l,T) + «^l,±(5r±,± - ffr:p,T)' ^tC. 
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Again by ()6.6p - ()6.8p . the initial conditions ()6.25|) arc uniquely determined by 
the initial divisor Vp^^no.to r) l|6-3p and by the data in (|6.5p . 

Being autonomous with polynomial right-hand sides, there exists a Tp > 0, such 
that the first-order initial value problem (|6.15p , (|6.23p , (|6.24p with initial conditions 
(|6.25p has a unique solution 

h,± = fL±{no,tr), £ = 0, . . . ,p± - 1, 

9L±'^ 9i,±{no,tr), e^O,...,p±, (6.30) 
hi^± = hi^±{no,tr_), ^ = 0, . . . ,p± - 1, 

for aU tr_ e {to,r - To, tf),r + To) 

(cf., e.g., [m Sect. III. 10]). Given the solution (|6.30p . we proceed as in Section|4] 
and introduce the following quantities (where tr G (io,r — lo, to,r + Tq)): 

, 2p+l . 1/2 9(no,tr) 

a+(no,t.) = a(no,^)( n -^™) IT Aifc(«o, (6.31) 

^ m=0 ' k=l 
P- P+-1 

= -co,+a+ino,tr)z-P- J] - Mfc(r^o, ^ )'"=^"'"*-\ (6.32) 

fe=i 

Gp{z,no,tr) = 7;(^JZ~1~^ ~ ^ ^ Fp{z,no,tr) (6.33) 

-P=(C:'j)=Afc("0,*r) 



2\a{no,tr) a+{nQ,tr)J - 



—Fp{z,no,tr) 



2a+(no,t.) ' ^' ^ ^ fJMno,U)~ £-1)1 



C=Mfc(no,tT:) 



In particular, this leads to the divisor 

2?A(no,tJ e SymP(/Cp) (6.34) 
and the sign of y in (|6.32p is chosen as usual by 

flk{no,tr) = (/ifc(rio,tr), {2/co, + )Hj{no,tr)'''Gp{Hk{no,tr),no,tr)), 

fc = 1, . . .,q(no,tr), 

and 

jl{no,tr) = {/il(no,tr), . . . ,/ll(no,tr), . . • ,Aig(«o,tr)("0,tr), . . . , Aig(„o (^0 , ) } 

" » ' . ' 

pi(na,t£j times Pg(7io ,tr) ("0>tr) times 

(6.36) 

with 

IJ^k{na,tr) ^ tJ.k'{no,tr) for A: 7^ fc', fc. A;' = 1, . . . , g(no, tr), (6.37) 

and 

Pk{no,tr) eN, k^l,...,q{no,tr), ^ pk{no,tr) = p. (6.38) 

fc=i 
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By construction (cf. ()6.35p ). the divisor P/i(„o.t,.) is nonspecial for all tr G (to,r — 
To, to^r + To). 

In exactly the same manner as in (|4.19p - (|4.2ip one then infers that Fp{- ,no, tr) 
divides Rp — Gp (since tr is just a fixed parameter). 

As in Section m the assumption that the Laurent polynomial Fp{- ,no — l,tr) is 
of full order is implied by the hypothesis that 

%no,tj e Mo for all tr e {to^r - To, to^r + To). (6.39) 

The explicit formula for /3(no, tr) then reads (for tr G (io.r — To, to,r + To)) 

a+{no,tr)(3{no,tr_) (6.40) 

^ 1 '^"^^^ (dP^("o.'J-i(C-^y(P))/dC^^("o-«^)-i)|^^^^,^^^,^^„^^^^) 

2 {pk{no,tr)-l)l 

q{no,tr) 

X n (Mfc ("0, tr)-Mfe' (no, t.))-P'= 

k' = l, k'^k 

-, / +/ , X 9("o,tr) 2p+l . 

1 / a^[no,tr) v-^ ,1 \ 

+ 7 ri~ + Pfc("-0,tr)/ifc(no,tr) - 77 > . (6.41) 

2ya{no,tr) fr[ - - ^,t^o ^ 

With (|6.2ip - (|6.4ip in place, we can now apply the stationary formalism as sum- 
marized in Theorem 14.41 subject to the additional hypothesis ()6.39p . for each fixed 
tr G (^o,r — To,to.r + To). This yields, in particular, the quantities 

Fp, Gp, Hp, a, (3, and /i, v_ for (ri, tr) el^x (to,r — To, to,r + To), (6.42) 

which are of the form (|6.32p - (|6.4ip . replacing the fixed np £ Z by an arbitrary 
n G Z. In addition, one has the following result. 

Lemma 6.2. Assume Hypothesis 16.11 and condition (|6.39p . Then the following 
relations are valid on <C x li x {to.r ~ To, to.r + To), 



Gp — 


FpHp 


— Rp, 


(6.43) 


z{Gp - Gp) + zpFp^ 


aHp 


= 0, 


(6.44) 


zl3Fp + aHp - Gp 


+ Gp 


= 0, 


(6.45) 


~Fp + zFp + a{Gp - 


f Gp) 


= 0, 


(6.46) 


z(3{Gp + Gp)- zHp 


+ Hp 


= 0, 


(6.47) 



and hence the stationary part, (|5.8p . of the algehro- geometric initial value problem 
holds, 

UVp - V+U = on C X Z X [to^r - To,to,r + Tq). (6.48) 
In particular. Lemmas I3.2fl3.5l apply. 



Lemma 16.21 now raises the following important consistency issue: On the one 
hand, one can solve the initial value problem (j6.57p . (j6.58p at n = no in some 
interval tr G (io.r ~ To, to.r + Tq), and then extend the quantities Fp,Gp, Hp to 
all C X Z X (to^r — To,to,r + To) using the stationary algorithm summarized in 
Theorem 14.41 as just recorded in Lemma 16.21 On the other hand, one can solve 
the initial value problem (|6.57p . (|6.58p a.t n — m, ni no, in some interval 
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tr G {to,r — Ti, to.r + Ti) with the initial condition obtained by applying the discrete 
algorithm to the quantities Fp, Gp, Hp starting at (jiq, tg.r) and ending at (rii, tg.r)- 
Consistency then requires that the two approaches yield the same result at n = ni 
for tr in some open neighborhood of io,r- 

Equivalently, and pictorially speaking, envisage a vertical i^-axis and a horizontal 
n-axis. Then, consistency demands that first solving the initial value problem 
(|6.57p . (|6.58p at n = no in some i^- interval around to.r and using the stationary 
algorithm to extend Fp,Gp, Hp horizontally to n = ni and the same t^-interval 
around io,r, or first applying the stationary algorithm starting at (no, io,r) to extend 
Fp, Gp, Hp horizontally to (ni, io,r) and then solving the initial value problem (|6.57p . 
(|6.58p at n = ni in some i^-interval around to,r should produce the same result at 
n = ni in a sufficiently small open tr interval around io,r- 

To settle this consistency issue, we will prove the following result. To this end 
we find it convenient to replace the initial value problem (|6.57p . (|6.58p by the 
original ir-dependent zero-curvature equation (|5.7p . Ut^ + UVr_ — V^U = on 

c X z X (to,.-ro,to,, + ro)- 

Lemma 6.3. Assume Hvvothesis \6.1\ and condition (j6.39p . Moreover, suppose that 
([QT l - dOOl) hold onCx {no} x {to^r - To,to,r + To). Then (jOT)) - ((5:49| hold on 
C X Z X (to^r — Tq, to.r + Tq) , that is, 

Fp,trXz, n, tr) = ~2iGp{z, n, tr)Fr{z, n, tr) 

+ i{Gr{z, n, tr) + Kr{z, U, tr))Fp{z, n, tr), 

Gp.triz, n, tr) = iFp(z, n, tr)Hr{z, n, tr) — iHp{z, n, tr)Fr{z, n, tr), 
Hp^tr-{z, n, tr) = 2iGp{z, n, tr)Hr{z, n, tr) 

— i(Gr{z, n, tr) + Kr{z, U, tr))Hp{z, U, tr), 
(Z, n,tr)eCxZx {to,r - To, to,r + Tq). 

Moreover, 

4>tA.P, n, tr) = iFriz, n, tr)<j)'^ (P, n, tr) 

- i(Gr{z, n, tr) + Kr_{z, n, tr_))(j){P, n, tr) + iHr{z, n, tr), 
cttjn, tr) = izFr{z, n, tr) 

+ ia(n, tr){Gr_{z, n, tr) + Kr {z, n, tr)) — iFr_{z, n, tr), 
Ptr {"n, tr) — —i(3{n, tr)(Gr {z, n, tr) + Kr{z, n, tr)) 
+ iHr{z, n, tr) — iz^^Hr {z, n, tr), 
(z, n, t,) e C X Z X {to.r - To, to,r + To). 

Proof. By LemmaEJwc have ({5:22]) . (lOl . (15:291) . ([OTIl - llOa)) . and ((OSll - fOTi) 

for (n, ) G Z X {to^r — To, to,r + To) at our disposal. 

Differentiating (I5.22p at n = no with respect to tr and inserting ()6.49p and ()6.50|) 
at n = no then yields (I6.52p at n = no. 

We note that the sequences fi^±, ge,±, hi^± satisfy the recursion relations (|2.6p - 
(|2.12p (since the homogeneous sequences satisfy these relations). Hence, to prove 



(6.49) 
(6.50) 

(6.51) 



(6.52) 
(6.53) 
(6.54) 



38 



F. GESZTESY, H. HOLDEN, J. MICHOR, AND G. TESCHL 



()6.53p and (|6.54|) at n — uq it remains to show 

. (6.55) 

= -iP{9r+,+ +9r--) - 'i{K+-l,+ - K--1-)- 

But this follows from ([09ll . (|63T|l at n = no (cf. (|6l8l) . ([QOp 

= - + «(/r_ - - /r+-l, + ), 

Inserting now (PTTTjl at £ = r_ - 1 and ([Ml) at £ = r+ - 1 then yields 

For the step n = no T 1 we differentiate (|4.36l) - (l4.4ip (which are equivalent 
to ([051) - ([071) ) and insert ([Ogi) - (p3T|) . ((?TT|) - ((ETa at n = no. For the case 
n > no we obtain and (if from (I6.49p . (|6.51|) at n = no as before using the 
other two signs in (|6.18p . (|6.29p . Iterating these arguments proves (|6.49p - (|6.54p for 
(z, n, i^) e C X Z X (io.r - To, io.r + To). □ 

We summarize Lemmas 16.21 and 16.31 next. 

Theorem 6.4. Assume Hypothesis 16.11 and condition (|6.39p . Moreover, suppose 
that 

h,± = fL±(nQ,tr), e = 0,...,p±-l, 

9i,± ^ 9L±{na,tr), e^O,...,p±, (6.56) 
hi^± = hi^±{no,tr), £^0,...,p±-l 

for all tr e {to,r - Jo, to,r + Jq), 

satisfy the autonomous first- order system of ordinary differential equations (for fixed 
n = no) 

fi,±,u = ^i,±{fk-,fk,+,9k-,9k,+), ^ = 0, . . . ,p± - 1, 

9£,±.t^ = Ge,±{fk-, fk,+ ,hk,~,hk,+)., e^O,...,p±, (6.57) 

hi,±.t^ = 'Ht±{gk-,9k,+, hk-,hk,+), i ^ 0, . . . ,p± - 1, 

with J-i,±, Ge,±, Ti.i^± given by (I6.17p . (|6.26p . (j6.27p . and with initial conditions 

fi,±{no,to,r), £ = 0, . . . ,p± - 1, 

9f.,±{'no,to^r), £^0,...,p±, (6.58) 
hi,±ino,to.r), ^ = 0, . . . ,p± - 1. 

Then Fp, Gp, and Hp as constructed in ()6.32p - ()6.42p on C x Z x (<o,r — To, to,r + Jo) 
satisfy the zero- curvature equations (|5.7p . (|5.8p . and (|5.50p on C x Z x (to.r — 
Jo, io,£ + Jq), 



Ut^ + UVr--Vf;U = 0, (6.59) 

c/yp-v;,+c/ = o, (6.6O) 

Vp,u-[Vr,Vp]^0 (6.61) 
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with U, Vp, and Vr given by (15. 9p . In particular, a, (3 satisfy (jS.ip and the algebro- 
geometric initial value problem (|5.2p . (j5.3[) on Z x (io.r — TQ,tQ,r + To)^ 

y-«/3t,+/3(.9^^,+ +5,._,-)-/ir_-i,-+/i.r+_i^+y (6.62) 

In addition, a, (3 are given by 

,2p+i . 1/2 g(",tr) 



2 



a+{n,tr)l3{n,tr) ^ -- }^ ( i ^ \ T\i 

X n (6.65) 

fe' = l, k'=/^k 

,,2p+l s l/2'?(".t) 9(n,tr) 

91 ( n^™) n MfcKir)"'"=^"'*-'+ 51 pfeK^rKKy 

m=0 ^ /c=l fe=l 

1 ^^"^^ \ 

- - ^ j , (z, n, ij.) e Z X (<o,j. - To, <o,r + To). 

m=0 ^ 

Moreover, Lemmas I3.2H3.5I and I5.2H5.4I apply. 

As in Lemma l4.3l we now show that also in the time-dependent case, most initial 
divisors are well-behaved in the sense that the corresponding divisor trajectory 
stays away from Poo± , Po,± for all {n, t^) S Z x M. 

Lemma 6.5. The set A4i of initial divisors T>fi,(na.tor) f'^'"' which Dfi^n.t^) o,''^d 
defined via (j5.58p and (|5.59[) . are admissible (i.e., do not contain Poo±, 
Po,±) <md hence are nonspecial for all {n,tr) € Z x M, forms a dense set of full 
measure in the set Sym''(/Cp) of nonnegative divisors of degree p. 



Proof. Let Alsing be as introduced in the proof of Lemma 14731 Then 



U [ag^iMsing) + trU^ 



^ U U (4Qo(^)+aQ„(Sym^'-i(/Cp)) + i.£f ) (6.66) 

-Pe{-Poo±,Po,±}tr:eR 

is of measure zero as well, since it is contained in the image of R X SymP-\jCp) 
which misses one real dimension in comparison to the 2p real dimensions of J{1Cp). 
But then 

U {^Qo (-^-ng) + nAp^^ _ (Poo+ ) + tjj^l^ 
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u( U (tiQo(-^si„g)+nAp„.jPooJ+t.C/r)+4p„.jPooj) (6.67) 

is also of measure zero. Applying to the complement of the set in (|6.67p then 
yields a set A^i of full measure in Sym'' (/Cp). In particular, A^i is necessarily dense 
in Sym^ (/Cp). □ 

Theorem 6.6. Let ^?/i(rio,to G M.\ he an initial divisor as in Lemma \6.5\ Then 
the sequences a, (3 constructed from /i(rio,io,r) as described in Theorem KAl satisfy 
Hvvothesis \5 .![ In particular, the solution a, (3 of the algebro- geometric initial value 
problem (pTM)) . (p^M)) is global in {n,tr) G Z x M. 

Proof. Starting with 'Djx[no,tar) ^ -^i: the procedure outlined in this section and 
summarized in Theorem [63] leads to 'Djx{n.t^) and T>{,i^n,t.r) for all (n, e Z x (io,r — 
7o,io,r + To) such that ((5?58)) and ((5759)) hold. But if a,/3 should blow up, then 
T^ii{n,tr) T^o{n,tr) must hit one of Foo± or Po,±, which is excluded by our choice 
of initial condition. □ 

We note, however, that in general (i.e., unless one is, e.g., in the special periodic 
case), 'Djj^(^n,u) S^t arbitrarily close to Poo±, Po,± since straight motions on 
the torus are generically dense (see e.g. T^, Sect. 51] or [551 Sects. 1.4, 1.5]) and 
hence no uniform bound (and no uniform bound away from zero) on the sequences 
a(n, tr)^ (3{n, tr) exists as (n, t^) varies in Z x M. In particular, these complex-valued 
algebro-geometric solutions of the Ablowitz-Ladik hierarchy initial value problem, 
in general, will not be quasi-periodic with respect to n or tr (cf. the usual definition 
of quasi-periodic functions, e.g., in [461 P- 31]). 



Appendix A. Hyperelliptic Curves in a Nutshell 

We provide a very brief summary of some of the fundamental properties and 
notations needed from the theory of hyperelliptic curves. More details can be 
found in some of the standard textbooks [50] , [21] , and [12] , as well as monographs 
dedicated to integrable systems such as [13l Ch. 2], [29| App. A, B], ^5V, App. A]. 

Fix p €N. The hyperelliptic curve ICp of genus p used in Sections [3]-[6| is defined 

by 

2p+l 

ICp-. Tp{z,y) = y"^ - R2p+2{z) ^0, i?2p+2(2) = J]^ (z - i^m), (A.l) 

m=0 

{Em}m=o,...,2p+i C C, Em ^ E„i' for ui ^ m' , m,m' = 0, . . . ,2p+ 1. (A. 2) 

The curve (jA.ip is compactificd by adding the points Poo+ and Poo_ , Poc+ Poo- , 
at infinity. One then introduces an appropriate set of p + 1 nonintcrsecting cuts Cj 
joining Em(j) and Em'{j) and denotes 

c= y Cj, qnCk^d), j^k. (A.3) 
je{i,...,p+i} 



Defining the cut plane 



(A.4) 
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and introducing the holomorphic function 

/2p+l \ 1/2 



'■p+1 \"'" 
n - E'u) (A.5) 

T1=0 / 



on n with an appropriate choice of the square root branch in (|A.5p . one considers 

Mp = {iz,aR2p+2{z)'/^)\zeC, ae{±l}}U{Poo+,Poo_} (A.6) 

by extending i?2p+2( • to C. The hyperelhptic curve ICp is then the set A4p with 
its natural complex structure obtained upon gluing the two sheets of A4p crosswise 
along the cuts. The set of branch points B{ICp) of ICp is given by 

B{ICp) = {(£;,„, 0)}„=o,...,2p+i (A.7) 

and finite points P on ICp are denoted by P = iz,y), where y{P) denotes the 
meromorphic function on ICp satisfying J-p{z, y) — if' — i?2p+2(^) = 0. Local coor- 
dinates near Pq = (zq, yo) e ICp\ {B{ICp) U {Poo+ , Poo_ }) are given by (p^ = z - zq, 
near Poo± by Cp^^ = l/-^, and near branch points (i?moiO) G B(ICp) by Qb^q.o) = 
{z — EmoY^"^ ■ The Riemann surface ICp defined in this manner has topological genus 
p. Moreover, we introduce the holomorphic sheet exchange map (involution) 

*:/Cp^/Cp, P=[z,y)^P* = {z,~y), P^^^P*^^^P^^. (A.8) 

One verifies that dz/y is a holomorphic differential on ICp with zeros of order 
p — 1 at Poo± and hence 

Vj = — - — , J = (A.9) 

form a basis for the space of holomorphic differentials on ICp. Introducing the 
invertible matrix C in C^, 



C — (Cj,fc)i,fc=i,...,p, Cj.k — rjj, 

1-1 



(A.IO) 

c(k) = {ci{k),. . . ,Cp{k)), Cj{k)=^C^J^, j, fc=l,...,p, 

the corresponding basis of normalized holomorphic differentials ujj, j = 1, . . . ,p, on 
ICp is given by 

^Cj{e)rii, / ujj = 6j^k, j,k = l,...,p. (A.II) 



/ .'-3 



Here {oj, 6j}j=i,...,p is a homology basis for ICp with intersection matrix of the cycles 
satisfying 

flj o 6fc = (5j,fc, o flfe = 0, 6j o 6fc = 0, j, /c = 1, . . . ,p. (A. 12) 

Associated with the homology basis {^j , &j}j=i,...,p we also recall the canonical 
dissection of ICp along its cycles yielding the simply connected interior ICp of the 
fundamental polygon dICp given by 

dICp = aibia^^b^^a2b2a2^b2^ ■ ■ ■ a~^b~^ . (A. 13) 

Let A4{ICp) and Ai^{ICp) denote the set of meromorphic functions (0-forms) and 
meromorphic differentials (1-forms) on ICp. Holomorphic differentials are also called 
Abelian differentials of the first kind. Abelian differentials of the second kind, 
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(^(2) g Ai^{ICp), arc characterized by the property that all their residues vanish. 
They will usually be normalized by demanding that all their a-pcriods vanish, that 
is, cj^^-' = 0, j = 1, . . . ,p. Any meromorphic differential cu'^^^ on ICp not of the 
first or second kind is said to be of the third kind. A differential of the third kind 
jj(3) g M^{}Cp) is usually normalized by the vanishing of its a-pcriods, that is, 
J^, w^^^ = 0, j = 1, . . . ,p. A normal differential of the third kind (j^p^^p^ associated 

with two points Pi, P2 & ICp, Pi P2, by definition, has simple poles at Pj with 
residues (—1)-'"'"^, j = 1,2 and vanishing a-periods. 
Next, define the matrix r = {Tj^e)j^e=i^,,,^p by 



Tj,e = / '^j^ j,i=l,...,p. (A.14) 

Jbe 

Then 

Im(r) > and tjj = T£j, j, £ = 1, . . . ,p. (A. 15) 

Associated with r one introduces the period lattice 

Lp = {zeCP\z = m + nT, m,neZP}. (A.16) 

Next, fix a base point Qo G ICp\ {Po,±, PcyD±}, denote by J{K,p) = C^/Lp the 
Jacobi variety of Kp, and define the Abel map -Aq^ by 

Aq^- ICp^ J{ICp), Aqo{P)=( [ ^i,---, [ (modLp), Pe/Cp. 

(A.17) 

Similarly, we introduce 

Oq^ : Div(/Cp) ^ J(/Cp), V ^ {V) = ^ V{P)Aq^ (P), (A.18) 

PeKp 

where Div(A^p) denotes the set of divisors on /Cp. Here P: /Cp ^ Z is called a 
divisor on /Cp if I?(P) 7^ for only finitely many P e /Cp. (In the main body of this 
paper we will choose Qo to be one of the branch points, i.e., Qo G '^(/Cp), and for 
simplicity we will always choose the same path of integration from Qq to P in all 
Abelian integrals.) 

In connection with divisors on /Cp we shall employ the following (additive) no- 
tation, 

1>QoQ = ^Qo + ^Q. Vq=Vq,+--- + Vq^ , (A.19) 
Q = {Qi,...,g„} e Sym"/Cp, Qo G /Cp, to e N, 

where for any Q e A^p, 

and Sym" K^p denotes the nth symmetric product of /Cp. In particular, Sym™/Cp 
can be identified with the set of nonnegative divisors < X> € Div(/Cp) of degree 

TO. 

For / G M{lCp) \ {0}, LO € M^{K.p) \ {0} the divisors of / and lo are denoted 
by (/) and (w), respectively. Two divisors V, £ £ Div(/Cp) are called equivalent. 
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denoted by P - f , if and only if P - f = (/) for some / e M{IC.p) \ {0}. The 
divisor class [D] of V is then given by [D] = {£ E Div(ACp) | £ ^ V}. We recall that 

deg((/)) = 0, deg((c.)) - 2{p - 1), / G M{ICp) \ {0}, u; e M\}Cp) \ {0}, (A.21) 

where the degree deg(X') of V is given by deg(2?) = J2peic ^i^)- customary 
to call (/) (respectively, {to)) a principal (respectively, canonical) divisor. 
Introducing the complex linear spaces 

C{V) = {feM{ICp)\f = QoY{f)>V}, r{V)=dimC{V), (A.22) 

C\V) = {Loe M\)Cp) I cj = or (cj) > V}, i{V) = AunC^{V), (A.23) 

with i{T>) the index of speciality of V, one infers that deg(2?), r(2?), and i{T>) only 
depend on the divisor class [D] of V. Moreover, we recall the following fundamental 
fact. 



Theorem A. 1. Let Pg e Sym^* /Cp, Q = {Qi, . . . , Qp}. Then, 

1 < i{Vq) = s 



(A.24) 



if and only if {Qi, . . . , Qp} contains s pairings of the type {P, P*}. ( This includes, 
of course, branch points for which P — P* .) One has s < p/2. 



Appendix B. Some Interpolation Formulas 

In this appendix we recall a useful interpolation formula which goes beyond the 
standard Lagrange interpolation formula for polynomials in the sense that the zeros 
of the interpolating polynomial need not be distinct. 



Lemma B.l ( 30J). Let p e N and Sp-i be a polynomial of degree p 
addition, let Fp be a monic polynomial of degree p of the form 



fe=i 



k=l 



Then, 



Sp^^iz) ^ Fp{z)j2 E 



k=l 1=0 ^^'^ ' 



dP- 



{z-cr 




C=p-k 



In particular, Sp^i is uniquely determined by prescribing the p values 

Sp-iink), S'p^j^iHk), ■ . ■ , Sl^j:j^^\pk), k^l,...,q, 

at the given points /ii. . . . , ^g. 

Conversely, prescribing the p complex numbers 



(0) (1) 



Apk-i) 



there exists a unique polynomial Tp_i of degree p - 



1, 



,q, 



q Pfc — 1 

Tp^,{z)=Fp{z)Y,Y. 



k=i e=o 



e\{pk-e-iy. 



1. In 



(B.l) 



(B.2) 



z e 



(B.3) 



(B.4) 



(B.5) 
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X 



k' = l,k'=^k 



, z e 



such that 

(B.6) 

We briefly mention two special cases of (jB.2p . First, assume the generic case 
where all zeros of Fp are distinct, that is, 

q=P. Pk^^, /Xfc 7^ /Xfc, for fc ^ fc', /c, fc' = 1, . . . ,p. (B.7) 

In this case (jB.2l) reduces to the classical Lagrange interpolation formula 

^^-[^^'-^ zee. (B.8) 

Second, we consider the other extreme case where all zeros of Fp coincide, that is, 
q=l, pi=p, Fpiz) = {z-fiir, zeC. (B.9) 
In this case (jB.2p reduces of course to the Taylor expansion of Sp-i around z = /ii, 

Sp-i{z) = J2 '-'^^^" {z-t^iY, zee. (B.IO) 

Appendix C. Asymptotic Spectral Parameter Expansions 

In this appendix we consider asymptotic spectral parameter expansions of Fp/y, 
Gp/y, and Hp/y, the resulting recursion relations for the homogeneous coefficients 

fi, (ji, and /if, their connection with the nonhomogeneous coefficients /f, gi, and 
hi, and the connection between ci^± and ci{E^^) (cf. (jC.3[) ). For detailed proofs of 
the material in this section we refer to [221, [32]. We will employ the notation 

E^'^{E^\...,Etp\,). (C.l) 

We start with the following elementary result (a consequence of the binomial 
expansion) assuming 77 £ C such that \'q\ < min{|i5o|~^, . . . , |£'2p+i|~^}: 

/2p+l \ 00 

\m=0 / fc=0 

where 

co{E) - 1, 

Jo,--- J2p + 1— U 

JoH [-i2p+i=fc 

The first few coefficients explicitly are given by 

^ 2p+l ^ 2p+l ^ 2p+l 

co{E) = l, ci{E)^^-Y,E^, C2{E)^- E„,,E„,,--Y,El, etc. 



m— mi ,7712—0 rn— 

mi <m2 



(C.4) 
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Next we turn to asymptotic expansions. We recall the convention y{P) = 
^^-P-i + o{CP) near Foo± (where ( = l/z) and y{P) = ±(co,-/co,+) + 0(C) 
near Po.± (where = z). 

Theorem C.l ([32]). Assume (|3.2p . s-ALp{a, f3) = 0, and suppose P = {z,y) S 
\ {-Poo+ , ^oo_ }• Then z^^Fp/y, z^^Gp/y, and z^^ Hp/y have the following con- 
vergent expansions as P Poo±, respectively, P Po,±, 



zP- Fp_{z) \TTZ,k+^'^\ P^Poo^, C = 1A, 

co,+ y \±ET=ok-<:'^ P^Po.±, C = z, 

zP_H^^(TET=ohL+C', P^Poo^, C=1A, 



z. 



(C.5) 
(C.6) 
(C.7) 



where — l/z {resp., C, — z) is the local coordinate near Poo± {resp., Po.±) and 
fi,±j ge,±j and h£,± are the homogeneous versions of the coefficients fi^±, gi,±, and 
hi^± introduced in (|2.16p - (|2.18p . 

Moreover, the Em-dependent summation constants Ci^±, £ — 0,...,p±, in Fp, 
Gp, and Hp are given by 

Q.± = co,±c,(^±i), e = 0,...,p±. (C.8) 
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